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Abstract 



In this thesis, we consider two different problems relevant to general relativity. Over 
the last few years, opinions on physically relevant singularities occurring in FRW 
cosmologies have considerably changed. We present an extensive catalogue of such 
cosmological milestones using generalized power series both at the kinematical and 
dynamical level. We define the notion of "scale factor singularity" and explore its re- 
lation to polynomial and differential curvature singularities. We also extract dynamical 
information using the Friedmann equations and derive necessary and sufficient condi- 
tions for the existence of cosmological milestones such as big bangs, big crunches, big 
rips, sudden singularities and extremality events. Specifically, we provide a complete 
characterization of cosmological milestones for which the dominant energy condition 
is satisfied. The second problem looks at one of the very small number of serious 
alternatives to the usual concept of an astrophysical black hole, that is, the gravastar 
model developed by Mazur and Mottola. By considering a generalized class of similar 
models with continuous pressure (no infinitesimally thin shells) and negative central 
pressure, we demonstrate that gravastars cannot be perfect fluid spheres: anisotropc 
pressures are unavoidable. We provide bounds on the necessary anisotropic pressure 
and show that these transverse stresses that support a gravastar permit a higher com- 
pactness than is given by the Buchdahl-Bondi bound for perfect fluid stars. We also 
comment on the qualitative features of the equation of state that such gravastar- like 
objects without any horizon must have. 
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Chapter 1 
Introduction 



This thesis is composed of three main chapters, and explores two separate problems 
in classical general relativity. 

• The first chapter assembles some basic notions in general relativity that will be 
exploited in the two following research chapters. 

• The second chapter is relevant to cosmology. It presents general definitions 
as a function of the scale factor of generic cosmological milestones (big bangs, 
sudden singularities...), and investigates their properties both at the kinematical 
and dynamical level. 

• The last chapter explores the properties of a gravastar like object — that would 
be an alternative to black holes without any horizon. 

The first problem deals with cosmology and singularities. In Friedmann-Robert son- 
Walker cosmologies, the physically relevant singularities had traditionally been thought 
to be restricted to the "big bang" or "big crunch". Even the most important "sin- 
gularity theorems" (e.g. Hawking and Penrose [Jj3 p 240]) do not specify which type 
of singularities are dealt with. However, in the last few years, additional different 
types of singularities have appeared in the literature, making the list of cosmological 
singularities more extensive — with big rips and sudden singularities added to the 
mix, as well as renewed interest in non-singular cosmological events such as bounces 
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and turnarounds. We define the notion of "scale factor singularity" and define generic 
cosmological milestones using generalized power series of the scale factor of the uni- 
verse a(t). Once this classification is in place, we can investigate properties of these 
cosmological milestones on kinematical and dynamical levels. 

• Can we find interesting relations between this notion of "scale factor singularity" , 
and polynomial and differential curvature singularities? 

• Introducing the Friedmann equations, is it possible to place constraints on whether 
or not the classical energy conditions are satisfied at the cosmological milestones? 

Since the classification of the cosmological milestones as generalized power series of 
the scale factor is extremely general, the corresponding results are to a high degree 
model-independent. We will demonstrate that there exists a particular class of singu- 
larities that do not imply a polynomial curvature singularity near the milestone and 
only an even smaller sub-class of them do not imply a differential curvature singularity. 
We will also determine, with a minimum of technical assumptions, the necessary and 
sufficient conditions for the occurrences of all these cosmological milestones. 

The second problem looks at an alternative model to black holes, the "gravastar 
model" , developed by Mazur and Mottola. Instead of having a star contracting and 
collapsing until matter arrives at a singularity at the centre, the gravastar model sug- 
gests that a gravitationally collapsing star would force spacetime itself to undergo a 
phase transition that would prevent further collapse. Thus the star would be trans- 
formed into a spherical quantum vacuum surrounded by a form of super-dense matter. 
In the Mazur and Mottola model, the gr avast ar-like objects are composed of three 
layers (interior de Sitter space, outer region consisting of a finite-thickness shell of stiff 
matter, exterior Schwarzschild vacuum geometry) and two infinitesimally-thin shells 
with surface densities a±, and surface tensions $±. 

• Is it possible to replace the thin shell completely with a continuous layer of finite 
thickness? 

• Is the pressure anisotropy (implicit in the Mazur-Mottola infinitesimally thin 
shell) a necessity for any gr avast ar-like objects? 
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• Is it really possible to build a gravastar-like objects using only perfect fluid with 
a continuous layer of finite thickness? Could a horizon or naked singularity form 
and why? 

Assuming pressure is continuous and differentiable, we will first analyze the resulting 
static geometry with the resulting isotropic TOV equation, and point out all the fatal 
problems arising with this situation. Secondly, we will turn to the resulting anisotropic 
TOV equation, and discuss when the anisotropy is necessary for the TOV equation 
to hold and what the resulting bounds are on transverse pressures in the anisotropic 
region. Finally, we will discuss what specific features the equation of state must have 
for an horizon-avoiding gravastar-like object to be sustained. Note that we derive 
these properties from an agnostic point of view as to the existence or non-existence of 
gravastars. 

To analyze these two problems, we have made use of the theory of classical general 
relativity, but we have not used a quantum field theory point of view. 

This thesis has three chapters but also four appendices. The first appendix ex- 
hibits all of the most important spacetime metrics used in this thesis, the last three 
appendices are papers published on work relating to this thesis, they were produced 
as a collaboration with Dr. Matt Visser, Petarpa Booserm, Tristan Faber and Silke 
Weinfurtner. At the time of writing, two papers ( "Effective refractive index tensor for 
weak-field gravity") [TT] and (Gravastars must have anisotropic pressures) [TH] have 
been published in Classical Quantum Gravity, and one paper (Necessary and sufficient 
conditions for big bangs, bounces, crunches, rips, sudden singularities, and extremality 
events) ^Hj has been accepted for publication in Classical Quantum Gravity. 



3 



i 



Chapter 2 



Some essential notions in General 
Relativity 

This chapter introduces some basic notions in general relativity and cosmology. This 
will serve as the starting point for the specific topics investigated in this thesis. The 
notions of general relativity presented in this chapter are non-exhaustive. 

In 1915, Einstein formulated the theory of general relativity, a theory of space, time, 
and gravitation. The new viewpoint it introduced on the nature of space and time often 
appears to be abstruse as it goes against some deeply ingrained, intuitive notions and 
requires some specific and sophisticated mathematic tools such as differential geometry. 

Later developments, and renewed interest in general relativity began in the 1960s, 
when the theory was related to other areas of physics and astronomy. The modern 
theory of gravitational collapse, singularities, and black holes was developed at this 
time. 

A deeper understanding of general relativity is another factor for renewed interest. 
Indeed, further understanding of the laws of nature are necessary to make progress 
toward the goal of developing a quantum theory of gravitation. Quantum gravity is 
the field of theoretical physics attempting to unify the theory of quantum mechanics, 
which three of the fundamental forces of nature, with general relativity, the theory 
of the fourth fundamental force: gravity. The theory of relativity in its own right 
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makes many remarkable statements concerning the structure of space and time and 
the structure of the gravitational field. 

One of the conceptual cores of general relativity is the Equivalence principle which 
comes in three different versions: 

1. Weak Equivalence Principle: "All bodies (subject to no force other than gravity) 
will follow the same paths given the same initial positions and velocities." This 
is also often called the Principle of Uniqueness of Free Fall. 

2. Einstein's (1907) Principle of Equivalence of gravitation and inertia: "All mo- 
tions in an external static homogeneous gravitational field are identical to those 
in no gravitational field if referred to a uniformly accelerated coordinate system." 

3. Strong Equivalence Principle: "At any event, always and everywhere, it is pos- 
sible to choose a local inertial frame (L.I.F) such that in a sufficiently small 
spacetime neighbourhood all (non-gravitational) laws of nature take on their fa- 
miliar forms appropriate to the absence of gravity, namely the laws of special 
relativity." 

All the three versions are correct, but it is the Strong Equivalence Principle (now 
often called the Einstein Equivalence Principle) which is now viewed as the most fun- 
damental. 

Another conceptual core of general relativity is the field equations. 

2.1 Notions of metric, geodesies, affine connexion, 
Killing vector field 

The bending of space and time nearby an object is called a gravitational field, and its 
potential can be described by a spacetime metric commonly written as: 



ds 2 = g a bdx a dx b . 
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A spacetime consists of a four- dimensional manifold M and a Lorentzian metric g a b 
defined everywhere on M. The metric has a Lorentzian signature (— , +, +, +). Indices 
a and b run from (usually denoted as the time direction) to 3. With the components 
of the metric g a b, we can construct a connexion known as the affine connexion or 
Christoffel symbol: 

r a 6c = l -g ad (g db ,c + g dc ,b + g b c, d ) (2.2) 

where g ad is the inverse metric. The connection indirectly represents the coordinate 
acceleration of a free falling particle in a gravitational field. 

A geodesic is the curved-space generalization of the notion of a "straight line" 
in Euclidean space. By definition, a straight line is the path of the shortest distance 
between two points. More technically, a straight line is a path which parallel transports 
its own tangent vector., i.e a curve whose tangent, T a , satisfies the equation: 

T a V a T b = (2.3) 

If we consider a parametric equation of the curve x a = x a (X), we obtain the geodesic 
equation: 

d 2 x a „„ dx b dx c 



If we take two events X\ and X 2 , define AX = X 2 — Xi, and define the quadratic 
form 

g(AX, AX) = (AXfg(AX), (2.5) 
we have the following definitions: 

• If g(AX, AX) < 0, the two events are timelike separated. It means that normal 
particles can successfully travel from one event to the other. 

• If g(AX, AX) = 0, the two events are lightlike separated. It means that one has 
to travel exactly at the speed of light to travel from one event to the other. 

• If g(AX, AX) > 0, the two events are spacelike separated. It means that one 
would have to travel faster than light to get from one event to the other. 
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In general relativity, timelike geodesies represent spacetime motions of free falling 
particles and null geodesies represent the paths of light rays. 



Another important mathematical concept is the one of Killing vector field. Some 
propositions on Killing vector fields are typically used when investigating properties 
related to stars (e.g. static and stationary spacetimes). 

A Killing vector field is a vector field on a Riemannian manifold that preserves 
the metric. Killing fields are the infinitesimal generators of isometries; that is, flows 
generated by Killing fields are continuous isometries of the manifold. Specifically, a 
vector field £ a is a Killing field if the Lie derivative £ with respect to £ a of the metric 
g a b vanishes: 



Equation ()2.7|) is also known as Killing's equation, thus, the necessary and sufficient 
condition that £ a be a Killing vector field is that it satisfy the Killing's equation. 

If a spacetime has a Killing vector, then we know we can find a coordinate system 
in which the metric is independent of one of the coordinates. By far the most useful 
fact about Killing vectors is the following lemma: 

Lemma 1. Let £ a be a Killing vector field and let A be an affinely parametrized geodesic 
with tangent u a . Then ^ a u a is constant along A. 

Lemma ^ can be interpreted as saying that every Killing vector field gives rise to a 
conserved quantity for particles and light rays. This conserved quantity enables one to 
determine the gravitational redshift in stationary spacetimes and is extremely useful 
for integrating the geodesic equation when symmetries are present. 



According to Riemann's theory of curved manifolds, the geometry of space-time is 
completely described by the metric tensor g a b, which has 10 algebraically independent 



A g a b = 0. 

In terms of the covariant derivative, this is 




V a & + V fe £a = 0. 



(2.7) 



2.2 Tensors 
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components at each event. 



2.2.1 Riemann tensor 

The Riemann tensor, and its derived tensors (Ricci, Weyl), are the only tensors that 
can be constructed from the metric tensor and its first and second derivatives: 

R a bcd = d c T a bd — ddT a bc + r e M r a ec — T e bc T a ed (2.8) 

The Riemann curvature tensor has 20 algebraically independent components at 
each event. The components of the Riemann tensor identically satisfy a differential 
equation (the Bianchi identity ()2.9|) ): 

Va-Rrfefcc + V ' c Rdeab + Vfe-Rrf eca = (2.9) 

The Riemann curvature tensor R a bcd is used to define other important tensors. 



2.2.2 Ricci tensor, Ricci scalar, Weyl tensor 

The Ricci tensor R ab f|2.1(J|) . the Ricci scalar R (I2.11|) and the Weyl tensor C abc d 
for n ^ 3 dimensions (J2.12)) . are defined by contractions, in a manner analogous to 
decomposing a matrix into trace and tracefree parts. 

Rab = g^Rdacb = R C acb (2-10) 



R = g ab R ab (2.11) 

Cabcd = Rabcd + ~ (fiW-Rcb + 9bcRda ~ 9a C Rdb ~ 9bdRca) 

n — 2 

+ 7 7T7 n\ {9ac9bd - gadgcb) R (2. 12) 

[n — l)(n — 2) 

If the Riemann tensor vanishes on a neighborhood of space-time, this neighbor- 
hood is locally isometric to Minkowski space-time, it is locally flat. Otherwise, if the 
Weyl tensor vanishes on a neighborhood of space-time, the neighborhood is locally 
conformally equivalent to Minkowski space-time. Thus, the Riemann, Ricci, and Weyl 
tensors all have geometric meaning independent of any physical interpretation. 
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2.2.3 Einstein tensor 



Finally, the Ricci tensor and the Ricci scalar can be used to define the Einstein tensor 
(jZH): 

Gab = Rab — ~9ab R (2-13) 

2.2.4 Energy- momentum tensor 

Now let's introduce the energy- momentum tensor T ab (also called the stress-energy 
tensor), this is a symmetric T 2 tensor which formulates the energy- like aspects of a 
system: energy density, pressure, stress, and so on. 
The stress-energy tensor is defined by: 



p 






IF 



(2.14) 



where p is the energy density (at rest in this coordinate system), F is the energy flux 
(generalization of the Poynting vector), and II is the stress tensor (pressure, shear 
etc.). 

The three most important energy-momentum tensors in general relativity are, 
namely, matter or dust, perfect fluid and the electromagnetic field. 



2.3 Einstein's field equations 

We can now write the Einstein's field equations (|2.15jl . which govern how the metric 
responds to energy and momentum. 

Gab = — 5— Tab (2-15) 

cr 

Here G is Newton's constant of universal gravitation, G = 6.67257x 10~ n m 3 kg" 1 s~ 2 . 
In other words, Einstein's field equation states that 



Rab — ~R 9ab = — 5— T ab , (2.16) 
2 cr 
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or equivalently 

R ab = ^(r ab -±g ab Ty (2.17) 

where T is the trace of the stress energy tensor. Note that often, the gravitational 
constant G and the speed of light c are set equal to one, this is called using "geometrized 
units" . 

Thus, the Ricci curvature is directly coupled to the immediate presence of matter 
at a given event. If there is no mass-energy at a given event, the Ricci tensor vanishes. 
If it were not for the Weyl tensor, this would mean that matter here could not have a 
gravitational influence on distant matter separated by a void. Thus, the Weyl tensor 
represents that part of space-time curvature which can propagate across and curve up 
a void. 

There are two particularly physically important solutions of Einstein's field equa- 
tion that represent two extremes of curvature. The Kerr vacuum solution, which 
models space-time outside a rotating body such as a star, has zero Ricci curvature 
but nonzero Weyl curvature at each event. The FRW spacetimes, which model the 
universe on a very large scale, has zero Weyl curvature but nonzero Ricci curvature at 
each event. 



2.4 Weak field theory 

When the gravitational field is "weak" , general relativity can be reduced to a linearized 
theory. The weakness of the gravitational field allows us to decompose the metric into 
the flat Minkowski spacetime (J2.19j) of special relativity plus a small perturbation, 



9ab — Vab 



h 



ab 



<1, 



:2.i8i 



where 



Vab 



-10 
10 
10 
1 



[2.19) 
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The assumption that h ab is small allows us to ignore anything that is higher than 
first order in this quantity, from which we immediately obtain: 



g ab = n ah + h ah + 0{h 2 ); h ab = rj ac ri bd h 



cd 



(2.20) 



The Christoffel symbol can now be rewritten as : 

1 



The Riemann tensor now becomes: 

1 



Rabcd — ~ 7^ (h ac ,bd + hbd,ac ~ ^ad,bc ~ hb c ,ad) + 0(Jl 



The Ricci tensor is: 



R ab = \ {K M + ht AC - v 2 h ab - h ab ) + 0(h 2 ). 



(2.21) 



(2.22) 



(2.23) 



By making coordinate changes, we can adopt the Einstein gauge (also called de Donder 
gauge, or Hilbert gauge, or Fock gauge) which provides the simple result, 



h c b - -hS c 
2 



0. 



and therefore, in this gauge, we have the simple results, 

1. 



R, 



fib 



-y 2 h ab + o(h 2 



(2.24) 



(2.25) 



and 



G 



= -^V 2 (hab - ^hVab ) 



Thus, using the weak field limit, we can write 



V 2 h ab = -16vrG T ab - -T Vab + 0(h 



(2.26) 



(2.27) 



which gives a lot more information. 
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2.5 Strong field theory 



Let's move from the domain of the weak field limit to solutions of the full nonlinear 
Einstein's equations. 

2.5.1 Black holes, event horizon 

Horizons are typically associated with strong gravitational fields (it is still possible to 
have a strong gravitational field without forming a horizon). By definition an event 
horizon is the boundary of the region from which future-directed null curves cannot 
escape to infinity. In other words, event horizons are "one-way membranes" permitting 
the passage of light and matter in only one direction and at which time slows to a 
stop. 

Some stars (whose masses are of the order of the sun's mass), can reach a final 
equilibrium state, e.g. white dwarf or a neutron star. However, for stars with much 
larger masses, such an equilibrium state is not possible. In this case, stars will contract 
to such an extent that the gravitational effects will overcome the internal pressure 
and stresses which will not oppose the contraction. The theory of general relativity 
predicts that such spherically symmetric stars will necessarily contract until all matter 
contained in the star arrives at a singularity at the centre of the symmetry. Those 
specific stars are surrounded by an event horizon and are referred to as black holes. 

2.5.2 Schwarzschild's solution to the Einstein equations 

One of the most important solution is the one discovered by Schwarzschild, which 
describes spherically symmetric vacuum spacetimes. Since it is in vacuum, Einstein's 
equations become: 



Thus, to obtain a solution, we need to find all four- dimensional Lorentz signature met- 
rics whose Ricci tensor vanishes and which are static and possess spherical symmetry. 
But first, we need to define more precisely the meaning of the terms "stationary", 
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"static" and "spherically symmetric" and to choose a convenient coordinate system 
for analyzing this class of spacetimes. 

• A stationary spacetime is a spacetime that possess a timelike Killing vector field 1 , 
£ a , (more precisely, a Killing vector field that is timelike near spatial infinity). 

• A spacetime is static if it is stationary and, if in addition, there exists a spacelike 
hypersurface £ orthogonal to the orbits of the isometry. It is equivalent to the 
requirement that the timelike Killing vector field £ a satisfy: 

£[aV fe £ c] = 0. (2.29) 

This implies that the spacetime metric ds 2 can be chosen to be invariant under 
a time revearsal about the origin of time, e.g. all cross terms dx t dx a V a / i 
vanish in the chosen coordinate system with arbitrary {x a } on E. Thus, in these 
coordinates, the metric components are of the form: 

3 

d s 2 = -V 2 (x\x 2 , x 3 )dt 2 + 9ab(x\x 2 , x 3 )dx a dx\ (2.30) 

a, 6=1 

where V 2 = — £ a £ a - The absence of dtdx a cross terms expresses the orthogonality 
of £ a with S. From the explicit form of a static metric ()2.30|) it is clear that the 
diffeomorphism defined by t — > — t is an isometry. And hence, in addition to 
the time translation symmetry, t — > t + constant, possessed by all stationary 
spacetimes, the static spacetimes also posses a time reflection symmetry. 
The static spherically symmetric metric will describe non-rotating stars or black 
holes, while rotating systems (which keep rotating in the same way at all times) 
will be described by stationary metrics. 

• A spherically symmetric solution means that there exists a privileged point, e.g. 
the origin, such that the system is invariant under spatial rotations about the 
origin. 

1 See section I^TI on Killing vector fields. 
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If the spacetime is stationary and spherically symmetric, it automatically implies that 
it is static. An example of a stationary but non-static case is the Kerr black hole, as 
the rotation rate is not changing but the surfaces orthogonal to the timelike Killing 
field are not spacelike. For stars we will always be looking at stationary cases (and 
therefore static by spherical symmetry). 

The Schwarzschild solution, describing static spherical symmetric vacuum space- 
times, is in standard coordinates (t,r,9,(p): 

ds 2 = - (l - 2 ^^j dt 2 + ^— \am dr 2 + rW (2.31) 

where dQ 2 = d6 2 + sin 2 9d(f> 2 . This is true for any spherically symmetric vacuum 
solution to Einstein's equations. M is a parameter that can be interpreted as the 
conventional Newtonian mass that would be measured by studying orbits at large 
distances from the gravitating source. 

As M — > we recover the expected Minkowski space. The metric also has 
the property of asymptotic flatness, as r — ► oo, the metric becomes progressively 
Minkowskian. 

Theorem 1. (Birkhoff's theorem) In a spherically symmetric spacetime the only so- 
lution to the Einstein equations (without cosmological constant) in vacuum is a piece 
of the Schwarzschild solution which is more properly called the "Schwarzschild exte- 
rior solution". Furthermore, a spherically symmetric vacuum solution in the exterior 
region is necessarily static. 

There is a generalization of Theorem [I] with a cosmological constant included which 
leads to the Kottler solution (Schwarzschild-de Sitter spacetime). It is interesting to 
note that the result is a static metric but that the source was not specified except that 
it be spherically symmetric. Specifically, the source itself does not need to be static, 
it could be a collapsing star, as long as the collapse were symmetric. Practically, the 
metric of the Schwarzschild exterior solution corresponds to (|2.31j) restricted to the 
region r > 2GM. 
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From the form of the Schwarzschild metric ([2.31)1 . we can see that the metric coef- 
ficients become infinite at r = and r = 2GM . The metric coefficients are coordinate- 
dependent quantities, but it is certainly possible to have a coordinate singularity which 
results from a breakdown of a specific coordinate system rather than the underlying 
manifold. The metric appears to have a singularity at r = 2GM, but there is actu- 
ally no physical singularity at that point, the coordinate system is just breaking down 
there. The value r = 2m is known as the Schwarzschild radius 2 , this is a removable 
coordinate singularity as indicated by the Riemann tensor scalar invariant 

Ra bcd R abcd = 48m 2 r" 6 , (2.32) 

which is finite at r = 2m. 

However there is an interesting physics phenomenon associated with this surface: 
it is the event horizon of the Schwarzschild black hole. 

2.5.3 Singularities 

The term singularity is used in many ways and in many theorems without being 
specifically, physically and fully explained. In this section, we will discuss some main 
definitions of what a singularity stand for and what it involves qualitatively. 

In a manifold M, there can be two points which are not connected by any causal 
curve (timelike or lightlike) (in a Euclidean signature metric this is impossible, but 
not in a Lorentzian spacetime). In such a case, it is said that the geodesic runs into 
a singularity, which can be thought of as the edge of the manifold. Manifolds which 
have such singularities are known as geodesically incomplete. In fact the "singularity 
theorems" of Hawking and Penrose ([7^1 p 240]) state that, for reasonable matter con- 
tent (no negative energies), spacetimes in general relativity are almost guaranteed to 
be geodesically incomplete. 

The Big bang or the initial singularity is defined by the fact that the scale factor 
of the universe a(t) — > in a finite time. This particular singularity results from a 
2 The value m is actually interpreted as the rclativistic mass with m = GM/c 2 in metres. 
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homogeneous contraction of space down to "zero size" , but does not represent an ex- 
plosion of matter concentrated at a point of preexisting nonsingular spacetime. A big 
rip is a singularity for which the scale factor of the universe a(t) — > oo in a finite time 3 . 

What kind of coordinate-independent feature is a warning that something goes 
wrong in the geometry, e.g a curvature singularity occurs? One simple criterion is 
when the curvature becomes infinite. The curvature is measured by the Riemann ten- 
sor from which can be constructed various scalar quantities, such as the Ricci scalar R, 
or higher order scalars R ab R a b, R abcd R a t, c d, etc... If any of those scalars go to infinity 
when approaching some point x, there is a curvature singularity at this point x. 

A naked singularity is any curvature singularity not surrounded by an event hori- 
zon. Since there is no event horizon, there is no obstruction to an observer travelling 
to the singularity and returning to report on what was observed. 

2.6 Energy conditions 

In classical general relativity, there are several types of energy conditions [73]: 

• the null energy condition (NEC); 

• the weak energy condition (WEC); 

• the strong energy condition (SEC); 

• the dominant energy condition (DEC). 

The energy conditions of general relativity permit one to deduce very powerful and 
general theorems about the behaviour of strong gravitational fields and cosmological 
geometries. In an orthonormal frame, the components of the stress energy tensor are 

3 For more definitions, see Chapter 
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given by: 





' P 
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rpab _ 





Pi 
















P2 







_ 








P3 . 



(2.33) 



The components of T are the energy density and the three principal pressures. 

2.6.1 Null Energy condition (NEC) 

For all future pointing null vectors k a , we ask that: 

T ab k a k b > (2.34) 
In terms of pressures and density, we have: 

Vi P + Pi ^ 0. 

Hawkings area theorem for black hole horizon relies on the NEC, and hence evap- 
oration of a black hole must violate the NEC. 

2.6.2 Weak Energy condition (WEC) 

Sometimes it is useful to think about Einstein's equations without specifying the theory 
of matter from which T ab is derived. This leaves us with a great deal of arbitrariness, in 
the absence of some constraints on T ab , any metric can satisfy the Einstein equations. 
The real concern is the existence of solutions to Einstein's equations with "realistic" 
sources of energy and momentum. The most common property that is demanded of 
T ab is that it represent positive energy densities — no negative masses are allowed. In 
a locally inertial frame this requirement can be stated as p = Too > 0. To turn this 
into a coordinate- independent statement, we ask that: 

T ab V a V b ^ V timelike vector V 

In terms of pressures and density, we have: 

p^O and Vzp + pj^O. 
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Any timelike vector can be a tangent to an observers world line. The WEC condition 
states that the energy density measured by any timelike observer is non-negative. It 
seems like a fairly reasonable requirement, and many of the important theorems about 
solutions to general relativity (such as the singularity theorems of Hawking and Penrose 
([ZOl P 240]) ) rely on this condition or something very close to it. Unfortunately it is 
not set in stone; indeed, it is straightforward to invent otherwise respectable classical 
field theories which violate the WEC, and almost impossible to invent a quantum field 
theory which obeys it. Nevertheless, it is legitimate to assume that the WEC holds in 
all but the most extreme conditions. 

2.6.3 Strong Energy Condition (SEC) 

For any timelike vectors V a , we ask that: 

(r ab - \goA v a v b > o 

where T is the trace of the stress-energy tensor: T = T ab g ab . 
In terms of pressures and density, we have: 

T = -p + ^Pt 

i 

Vz p + Pi ^ and p + p.- L ^ 0. 

i 

Note that the SEC implies the NEC, it does not imply the WEC. For example, matter 
with a negative energy density but sufficiently high pressures could satisfy the SEC 
but would violate the WEC. 

The Penrose-Hawking singularity theorem relevant to the cosmological singularity 
uses the SEC. 

2.6.4 Dominant Energy Condition (DEC) 

For any timelike vectors V a , we ask that: 

T ab V a V b ^ and that T ab V b is a future directed non-spacelike vector. 
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The DEC assumes that the WEC holds, and that for all future directed timelike 
vectors V a that T a bV b is a future directed non-spacelike vector. This ensures that the 
net energy flow does not exceed the speed of light. The dominant energy condition 
implies the weak energy condition and also the null energy condition, but does not 
necessarily imply the strong energy condition. 

In terms of pressures and density, we have: 

p ^ and V i — p ^ Pi ^ p. 

The dominant energy condition can be interpreted as saying that the speed of energy 
flow of matter is always less than the speed of light. 



2.6.5 Comments 

Note that the null energy condition implies the weak energy condition, but otherwise 
the NEC, the WEC and the SEC are mathematically independent assumptions. In 
particular, the SEC does not imply the WEC. It is stronger only in the sense that it 
appears to be a stronger physical requirement to assume equation ()2.6.3|) rather than 
equation (13X21) . Violating the NEC implies violating the DEC, SEC and WEC as well. 

The energy conditions are looking a lot less secure than they once seemed: 

• There are quantum effects that violate all of the energy conditions. 

• There are even relatively benign looking classical systems that violate all the 
energy conditions |74| . 

Hawkings area theorem for black hole horizon relies on the NEC, and hence evap- 
oration of a black hole must violate the NEC. 
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2.7 Cosmology 



2.7.1 Introduction 

Cosmology is the study of the dynamical structure of the universe considered as a 
whole. Contemporary cosmological models are based on the idea that the universe 
is, on average, the same overall. This is based on a very simple principle, called "the 
cosmological principle" , which is a generalization of the Copernican principle: 

The cosmological principle: at each epoch, the universe presents the same aspect 
from every point, except for local irregularities [HK] . 

When averaged over sufficiently large volumes the universe and the matter in the 
universe should be isotropic and homogeneous. 

• Isotropy states that space looks the same no matter what direction one looks at 
(direction independent). 

• Homogeneity is the statement that the metric is the same throughout the space 
(position independent). 

Astronomical observations reveal that the universe is homogeneous and isotropic 
on the largest scales. Traditionally this homogeneity has been assumed up to "small" 
fluctuations that are large enough to include clusters of galaxies. The scale at which 
homogeneity sets in is still not completely certain. Voids with diameters of order 10 8 
light years are ubiquitous, forming at least 40% of the volume of the universe [IB], [4*§] . 
and are typically surrounded by bubble walls containing galaxy clusters. The largest 
feature observed - the Sloan Great Wall jlU] - is 1.47 x 10 9 light years long. We simply 
assume homogeneity for some suitable defined cell size. 

When looking at distant galaxies, they seem to be receding from our galaxy. It ap- 
pears that the universe is not static, but changing with time. Thus most cosmological 
models are built on the fact that the universe is homogeneous and isotropic in space, 
but not in time. Observationally, the universe today is significantly different from the 
universe of 10 10 years ago, and radically different from the universe of 1.5 x 10 10 years 
ago. 
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2.7.2 "Cosmography" 

Simply by using the assumptions of isotropy and homogeneity, a cosmological model 
can be derived without yet using Einstein equations. This homogeneous and isotropic 
cosmological model is called the Friedmann- Robertson- Walker geometry (|2.35|) . (FLRW 
Friedmann-Lemaitre-Robertson- Walker geometry) and is given by: 

ds 2 = -dt 2 + a(t) 2 | - ^ 2 + r 2 [d6 2 + sin 2 6 dcf 2 ] j (2.35) 

where, a(t) is the scale factor of the universe. There are only three values of interest 
for the parameter k: 

• k = — 1, this corresponds to a negative curvature (for the hyperboloid) 

• k = 0, this corresponds to no curvature (flat space) 

• k = +1, this corresponds to a positive curvature (for the 3-sphere) 

Therefore, the assumptions of homogeneity and isotropy alone have determined the 
spacetime metric up to three discrete possibilities of spatial geometry k and the arbi- 
trary positive function of the scale factor a(t). 

Observational evidence strongly suggests that our universe (or part of our universe 
within our causal past), is well described by a Friedmann-Robertson- Walker model, 
and indeed a k = model, at least as far as the decoupling time of matter and 
radiation. 

2.7.3 "Cosmodynamics" 

Now, by substituting the spacetime metric (J2.35J) into Einstein's equations ()2.15|) . 
some predictions for the dynamical evolution of the system can be obtained. But first, 
we need to describe the matter content of the universe in terms of the stress-energy 
tensor. Using the assumptions of isotropy and homogeneity, the stress-energy tensor 
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of matter in the present universe is approximated in an orthonormal frame by: 





' P 
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rpab _ 
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V 
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p . 



(2.36) 



Here p and p are the average density and pressure due to the galaxies, stars, clouds of 
dusts etc... 

Applying the Einstein equations imply the two following Friedmann equations: 



a 2 



a- 



+ 



a" 



a 2 k a 

- + - + 2 - 

a z a z a 



8nG N p = 
8nG N p = 

And consequently, equation ()2.37j) and equation ()2.38j) imply: 

8ttG n [p + 3p] = -6-. 



(2.37) 
(2.38) 



(2.39) 



The Friedmann equations completely specify the evolution of the universe as a function 
of time. The difficulty remains to detemermine a suitable matter model for p and p, 
that is to make even more progress, it is necessary to choose an equation of state 
between p and p. 



2.7 A Cosmological parameters 

This section introduces some of the basic terminology associated with the cosmological 
parameters. 

The rate of expansion is characterized by the Hubble parameter: 



H 



a 
a 



(2.40) 



The Hubble parameter quantifies the "speed" with which the size of the universe 
is increasing. The value of the Hubble parameter at the present epoch is the Hubble 
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constant, H . There is currently a great deal of controversy about what its actual value 
is, currently the measurements 4 give: 

H = 100 h km/s/Mpc, (2.41) 

with a present day normalized Hubble expansion rate h, 

h = 0.7ll°;° 4 (2.42) 

"Mpc" stands for "megaparsec" , 1 Mpc = 3 x 10 24 cm. 

The universe is expanding, therefore we know that a > 0. From equation ()2.39j) we 
also know that a < when assuming that the pressure p and the density p are both 
positive. The universe must have been expanding at a faster and faster rate when 
going back in time. If we consider that the universe have always been expanding at 
the present rate, then at the time T = H~ l = a/a ago, the scale factor a would be 
null, a = 0. However, the expansion rate was actually faster, therefore, the time at 
which a = was even closer to the present. By assuming homogeneity and isotropy, 
general relativity makes the prediction that at a time less than H~ l ago, the universe 
was in a singular state. This singular point referred to as the big bang had an infinite 
density of matter and an infinite curvature of spacetime. 



The value of the Hubble parameter changes over time either increasing or decreasing 
depending on the sign of the deceleration parameter: 

ad 

q = -— • 2.43 

a z 

The deceleration parameter measures the rate of change of the rate of expansion. 
Different values, or ranges of values, of go correspond to different cosmological models. 
In principle, it should be possible to determine the value of go observationally. For 
example, for a set of identical supernovae within remote galaxies, the relationship 
between apparent brightness and redshift is dependent on the value of the deceleration 
parameter. Although measurements of this kind are notoriously difficult to make and 



1 See S. Eidelman et al. from the Particle Data Group J2H] for recent measurement values. 
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to interpret, recent observations tend to favor accelerating universe models. In [73], it 
was found that 

g = -0.551°;* (2.44) 
Another useful quantity is the energy density parameter, 

^ 8nG p . 

where the critical density (Hubble density) is 

Pcrit = (2.46) 

07TG 

This quantity (which will generally change with time) is called the "critical" density 
and current measurements 5 give: 

p crit = 2.775 x 10 11 h 2 M e Mpc~ 3 , (2.47) 

where M & is the solar mass and h is the present day normalized Hubble expansion 
rate. Using the Friedmann equation (|2.37|) . we can then write: 

- 1 = 77^- ( 2 - 48 ) 

The sign of k is therefore determined by whether the energy density parameter Q is 
greater than, equal to, or less than one. Indeed, 

P < Peru ^ Q < 1 «-> k = — 1 «-> open 

P = Pent *-> ^ = 1 ^ k = «-> flat (2.49) 

P > Pent ^ Q > 1 «-> fc = +1 «-> closed 

The density parameter, then, indicates which of the three Robert son- Walker geome- 
tries describes our universe. Determining it observationally is an area of intense inves- 
tigation, however, presently, it is thought to be 6 : 

Q = 1.02 ±0.02. (2.50) 



5 See S. Eidclman et al. from the Particle Data Group for recent measurement values. 
6 See S. Eidclman et al. from the Particle Data Group |29j for recent measurement values. 



25 



2.8 Conclusion 



This chapter recalls the main important features of general relativity that makes re- 
markable statements concerning the structure of space and time and the structure of 
the gravitational field (e.g. Einstein's equations). 

The energy conditions of general relativity permit one to deduce very powerful and 
general theorems about the behaviour of strong gravitational fields and cosmological 
geometries. 

A description of our universe has been outlined in section l2~71 Many important issues 
in cosmology remain to be solved, but it is already clear that general relativity has 
provided a successful picture of our universe. 
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Chapter 3 

Necessary and sufficient conditions 
for big bangs, bounces, crunches, 
rips, sudden singularities, and 
extremality events 

3.1 Introduction 



In section 12.71 we have discussed that when averaged over sufficiently large volumes 
the universe and the matter in the universe should be isotropic and homogeneous. 

As a consequence, the cosmological principle, which is ultimately a distillation of 
our knowledge of observational cosmology, leads one to consider cosmological space- 
times of the idealized FRW form [J7| HI EOJ Eg ■ 



where, a(t) is the scale factor of the universe and should be positive, and, where, the 
parameter k takes only three values: 

• k = — 1, negative curvature 




(3.1) 
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• k — 0, flat space 



• k — +1, positive curvature. 

The central question in cosmology is now the prediction of the history of the scale 
factor a(t) [13 HE]. 

As discussed in l2.7.3| and assuming we can apply the Einstein equations of classical 
general relativity to cosmology in the large, we get some key dynamical equations: the 
Friedmann equations f)3. 21 13.31 13.41 13. 5J) . 

Using units where 8ttGn = 1 and c = 1, we have: 



pit) 




a a 2 



k 



2 ■ 



(3.2) 



(3.3) 



a 



pit) + 3p(t) 




(3.4) 



and the related conservation equation 



p(t)a 3 + 3 [p(t) +p(t)]a 2 a 







(3.5) 



Theorem 2. Any two of equations HS. 6 A\S.'A \3.4\ ) imply the remaining one, and imply 
equation \3. 5)) . Equation \3. 5)) and any one of equations \3. 6 A \3.<A \3.J$ imply the two 
other equations but with a specific choice of integration constant. 



Lemma 2. Equations \3. gj) and \3. 3]) imply equations \3.4\l and LS.oV . 



Proof. 




(3.6) 



a 



(3d 2 + 3A; 



6da — 3d 2 — 3fc) 



(3.7) 
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Equation (|3.8|) is exactly equation (|3.4j) . Now, using equation (|3.2|) . we derive: 

3 

p(t) = — {2a 2 ixa - 2ad 3 - 2kaa) (3.9) 
ct 



p(t)a = Qada — 6d — Gka (3.10) 



We also have: 



3(p + p)a 2 a = 3d (-2aa + 2d + 2k) (3.11) 
= — 6add + 6d 3 + 6ka (3-12) 

Now adding equations ( 13.10)) and ()3.12|) : 

3 (p + p) « 2 d + p(t)a 3 = — 6add + 6d 3 + 6/cd + 6add — 6d 3 — 6ka (3.13) 

= (3.14) 

Therefore we obtain equation (|3.5J) from equations (|3.2J) and ( 13. 3|) . □ 
Lemma 3. Equations A3. SA) and \3.J^ imply equations A3. 3]) and A3. 51) . 
Proof. Using equation ()3.4j) . we can isolate p and replace p by equation ()3.2j) : 

> = -!-£■ < 3 - 15 > 

= -22-^-* (3.16) 

Thus, equations (|3.2j) and (|3.4J) imply equations (|3.3jl . 

And as proven for Lemma |2J, equations ()3.2j) and (|3.3|) imply equations (|3.4|) . □ 

Lemma 4. Equations A3. OA) and A3. 5\) imply equations A3. 3\) and \3.4\ j- 

Proof. Using equation ( |3.2|L we replace p and p in equation (|3.5|) to obtain an equation 
for the pressure p. 

p{t)a 3 + 3{p{t) + p{t))a 2 a = (3.17) 
==>- 3pa 2 a = — 9d 3 — 9ka — 6aaa + 6d 3 + 6kd (3.18) 

CL CL k , , 

=> p = -2 3.19 

a er a 1 

Thus, equations ()3.2j) and (|3.5jl imply equations (j3.3j) . 

And as proven for Lemma El , equations ()3.2|) and ()3.3|) imply equations (|3.4|) . □ 
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Lemma 5. Equations h3. 3\) and V3.5)) imply equations \3. ty) and up to a term 

proportional to an arbitrary integration constant K: p ~ ^. This integration constant 
has to be eliminated "by hand" 1 . 

Proof. Replacing the pressure p in equation (|3.5|h by its value given in equation ()3.3|) . 
we obtain the following differential equation: 

3pa 2 a + pa 3 = —3pa 2 a (3.20) 

ad a 3 ka , 
3pa + pa = 6— + 3^ + 3^- (3.21) 
a a z a z 

a ad a 3 ka , 

p + 3-p = 6-t + 3^ + 3^- 3.22 

a a z cr a 6 

By integrating this differential equation in p, we get: 

_ I { ex P (/ fdt) (6f + 3$ + 3§) } <ft + JT 



exp (/ f A) 

where K is an integration constant. Furthermore, 



and therefore, 



We finally have: 



(3.23) 



exp (^J ^d?j = exp (3 ln(a)) = a 3 , (3.24) 



-l>f*H^)}*= /W^l»* (3.25) 



y (6add + 3d 3 + 3kd) dt (3.26) 
3ad 2 + 3ka (3.27) 



fa 2 k K\ . 
The last term — - of equation ()3.28j) has to be eliminated "by hand" to obtain exactly 

3(2 

equation ()3.2|) . 

And as proven for Lemma EJ, equations ()3.2|) and ()3.3|) imply equations (|3.4j) . □ 



1 Physically, this extra term corresponds to "dust'' 
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Lemma 6. Equations and \3. 5)) imply equations \3. ty) and \3. 3\) with k arising 

as an integration constant and being set to k = — 1/0/ + 1 by suitably rescaling a(t). 

Proof. Using equations ()3.4|) and (|3.5J) . we get: 

p(t)a 3 + 3 ( p(t) + (-2--^-^X\a 2 a = (3.29) 

a ad , 

=> P + 2-p = 6— (3.30) 

a cr 

By integrating the differential equation, we obtain: 

[ \exp(f2^dt)6^\dt + K 

p= J 1 U "/ ■ "V , 3.31 

exp(/2fdfc) 

where K is an integration constant. Furthermore, 

exp (^J 2^dtj = exp (21n(a)) = a 2 , (3.32) 

J Qaddt + K 



and therefore, 



a 2 



/3f + K 



a 2 



3d 2 + K 



a 2 



(3.33) 
(3.34) 
(3.35) 



By suitably rescaling a(t), we can, without loss of generality set this integration 
constant to k = —1/0/ + 1. That is 

• if K = 0, then we do not need to rescale a(t). 

• if K 7^ 0, then we can define a(t) ncw = a(i)/ y/\K\. After this substitution, we 
get 



■ \f ' + K 



(3.36) 



dL„ K 



w new I- 11 - l^ncw 



3(-p + — (3.37) 



3(4^ + ^), (3-38) 



q1 q1 

new new 
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with the only possible values of k = 



-1/0/ + 1. 



Finally, we get the expected equation: 



) 



(3.39) 



And as proven for Lemma El, equations (|3.2jl and ()3.4j) imply equations ()3.3|) . □ 

The physically relevant singularities occurring in the Friedmann-Robertson- Walker 
cosmologies had traditionally been thought to be limited to the "big bang" , and pos- 
sibly a "big crunch". However, over the last few years, the zoo of cosmological sin- 
gularities considered in the literature has become considerably more extensive, with 
"big rips" and "sudden singularities" added to the mix, as well as renewed interest in 
non-singular cosmological events such as "bounces" and "turnarounds" 2 . 

Typically, those singularities are classified in relation to what happens to the scale 
factor a(t) of the universe at the time when the singularity occurs. 

• For a big bang and a big crunch, a(t) — ► 0. 

• For a big rip, a(t) — ► oo. 

• For a sudden singularity, the definition varies according to different authors: a(t) 
is finite at the time of the occurrence of the singularity (everyone agrees on that), 
however, the first derivative of the scale factor is infinite a(t) — ► oo, or it could 
be the second derivative a(t) — > oo, etc., depending on how "sudden" one is 
willing to describe a sudden singularity. 

• For bounces and turnarounds, a(t) is finite and reaches either a minimum or a 
maximum. 

• For points of inflexion, a(t) is finite and reaches a stationary point. 

2 See references |H1 [3 EI 113 CSI HHI EH ESI E2I E31 EH1 E2I EHI ESI Examples of sudden singularities 
have been given earlier, for example in references |55j . |56| . in which they are referred to as "Crack 
of Doom" singularities. 
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Problem: 

Is it possible to come up with generic definitions that would apply to a complete 
catalogue of such cosmological milestones, both at the kinematical and dynamical 
level? 

• Kinematics: 

Can we define a notion of "scale-factor singularity" and find interesting relations 
between this notion and curvature singularities? 

• Dynamics: 

Using the Friedmann equations (without assuming even the existence of any 
equation of state) is it possible to place constraints on whether or not the classical 
energy conditions are satisfied at the cosmological milestones? 

In this chapter, we use these considerations to derive necessary and sufficient con- 
ditions for the existence of cosmological milestones such as bangs, bounces, crunches, 
rips, sudden singularities, and extremality events. Since the classification is extremely 
general, the corresponding results are to a high degree model-independent: in par- 
ticular, we determine, with a minimum of technical assumptions, the necessary and 
sufficient conditions for the occurrences of all theses cosmological milestones. Note 
that we only use classical general relativity to derive those results. 

3.2 Cosmological milestones and kinematics. 
3.2.1 Definitions 

In this section, we develop precise definitions of a hopefully complete catalogue of 
cosmological milestones. 

It is clear that it would be really convenient to have some unspecified generic 
cosmological milestone defined in terms of the behaviour of the scale factor a(t), and 
which occurs at some finite time t & . 
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Generalized Frobenius series are commonly used when expanding solutions of dif- 
ferential equations around their singular points. Motivated by this property, we will 
assume that in the vicinity of the milestone the scale factor has a generalized power 
series expansion that generalizes the notions of Taylor series, meromorphic Laurent 
series, Frobenius series, and Liapunov expansions jHl] and are even more general than 
the generalized Frobenius series adopted in |78j . Indeed, in the present context, if the 
scale factor a(t) is representable by such a generalized power series, then by the Fried- 
mann equations both p(t) and p(t) are representable by such power series. Formal 
reversion of the power series then implies that the equation of state p(p), and thence 
also the function p(a) possess such generalized power series. Conversely, if p(a) is rep- 
resentable by such a generalized power series, then by the first Friedmann equation, 
a(t) has a power series of this type, which upon integration implies that a(t) itself 
possesses such a power series. Similarly, if the equation of state p(p) is representable 
by such a generalized power series then by integrating the conservation equation we 
have 



which will now also possess such a generalized power series. That an extension of the 
usual concept of a Frobenius series is likely to be useful is already clear from the analysis 
of [ZHj- We should also be clear concerning what type of object falls outside this class of 
generalized power series: First, essential singularities [effectively poles of infinite order, 
arising for example in functions such as exp(— considered in the neighbourhood 
of x = 0] lie outside this classification, and secondly, certain variations on the notion 
of Puiseux series [specifically, series containing (lnx) n , (lnlnx) n , (lnlnlnx)",...] also 
lie outside this classification. We are not aware of any situations in which these 
exceptional cases become physically relevant. 

Definition 1. Generic finite-time cosmological milestones: Suppose we have some 
unspecified generic cosmological milestone, that is defined in terms of the behaviour of 
the scale factor a(t), and which occurs at some finite time t Q . We will assume that 
in the vicinity of the milestone the scale factor has a (possibly one-sided) generalized 




(3.40) 
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power series expansion of the form 

a(t) = c |t - t p + Cl \t - top + c 2 \t - t p + c 3 |t - top + . . . (3.41) 

where the indicia! exponents rji are generically real (but are often non-integer) and 
without loss of generality are ordered in such a way that they satisfy 

7]o < r/i < 772 < 7/3 . . . (3.42) 

Finally we can also without loss of generality set 

c > 0. (3.43) 

There are no a priori constraints on the signs of the other Ci, though by definition 

From a physical point of view, this definition is really generic and can be applied 
to any type of cosmological milestone. This generalized power series expansion is suf- 
ficient to encompass all the models we are aware of in the literature, and as a matter 
of fact, the indicial exponents rji will be used to classify the type of cosmological mile- 
stone we are dealing with. For many of the calculations in this chapter, the first term 
in the expansion is dominant, but even for the most subtle of the explicit calculations 
below it will be sufficient to keep only the first three terms of the expansion: 

a(t) = co|t-t p+ci|t-top+c 2 |t-top...; Vo < Vi < V2; c > 0. (3.44) 

The lowest few of the indicial exponents are sufficient to determine the relationship 
between these cosmological milestones, the curvature singularities and even the energy 
conditions of classical general relativity. 

Note that this expansion fails if the cosmological milestone is pushed into the 
infinite past or infinite future. We shall return to this point later when we discuss the 
total age of the universe. 
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Big Bangs and big crunches 

Big bangs and big crunches are some of the most basic cosmological milestones en- 
countered in the literature. For these types of singularities, the scale factor a(t) — > 
at some finite time as one moves to the past or future. 

According to the Big Bang theory, the universe originated in an infinitely dense 
and physically paradoxical singularity. Space has expanded with the passage of time, 
objects being moved farther away from each other. In cosmology, the Big Bang theory 
is the prevailing theory about the early development and shape of the universe. The 
central idea is that the observation that galaxies appear to be receding from each other 
can be combined with the theory of general relativity to extrapolate the conditions of 
the universe back in time. This leads to the conclusion that as one goes back in time, 
the universe becomes increasingly hot and dense. However, the big bang represents 
the creation of the universe from a singular state, not explosion of matter into a 
pre-existing spacetime. It might be hoped that the perfect symmetry of the FRW 
universes was responsible for this singularity, but however, the singularity theorems 
predict (under relatively mild conditions) that any universe with p > and p ^ 
must have begun at a singularity. 

There are a number of consequences to this view. One consequence is that the 
universe now is very different than the universe in the past or in the future. The 
Big Bang theory predicts that at some point, the matter in the universe was hot and 
dense enough to prevent light from flowing freely in space. That this period of the 
universe would be observable in the form of cosmic background radiation (CBR) was 
first predicted in the 1940s, and the discovery of such radiation in the 1960s swung most 
scientific opinion against the Big Bang theory's chief rival, the steady state theory. 

In cosmology, the Big Crunch is a hypothesis that states the universe will stop 
expanding and start to collapse upon itself; a counterpart to the Big Bang. If the 
gravitational attraction of all the matter in the observable horizon is high enough, 
then it could stop the expansion of the universe, and then reverse it. The universe 
would then contract, in about the same time as the expansion took. Eventually, all 
matter and energy would be compressed back into a gravitational singularity. It is 
meaningless to ask what would happen after this, as time would stop in this singular- 
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ity as well. 



We shall define the scale factor near a big bang and a big crunch using the gener- 
alized power series expansion mentioned earlier. 

Definition 2. Let the time of the big bang (if one occurs) be denoted by and the 
time of the big crunch (if one occurs) be denoted by t®. We shall say that the bang or 
crunch behaves with indicial exponents (0 < 770 < 771 . . . ) if the scale factor possesses 
a generalized power series in the vicinity of the singularity. That is, if 

a(t) = c (t - + ci(t - i*)" 1 + . . . (3.45) 

or 

a(t) = c (t© - tf° + ci(t® - ty 1 + ... (3.46) 

respectively. Moreover, the series have been carefully constructed to make a(i*) = 
and a(t@>) = 0. 

Big rips 

A "big rip" is a singularity for which a(t) — * 00 at finite time. A big rip could occur in 
the future or in the past, however, the literature to date has solely considered future 
rips (as a past rip would be a most unusual and unexpected beginning to the history 
of the universe). 

The Big Rip is a cosmological hypothesis about the ultimate fate of the Universe. 
The key to this hypothesis is the amount of dark energy in the universe. If the universe 
contains enough dark energy, it could end with all matter being pulled apart. First 
the galaxies would be separated from each other, then gravity would be too weak to 
hold individual galaxies together. 

We shall write the scale factor near a big rip using the generalized power series 
expansion. 

Definition 3. Let the time of the rip, if it occurs, be denoted t%., then we define the 
indicial exponents of the rip (either future or past) to be (r]o < rji . . . ) if the scale factor 
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possesses a generalized power series in the vicinity of the rip: 

a(t) = cb|t* - t|"° + ci|t» - t\ m + ..., (3.47) 

with t]q < and Cq > 0. Moreover, the series has been carefully constructed to make 
a(t%.) = oo. 

Note the similarity to bangs and crunches, with the only difference being in the 
sign of the exponent rjo- 

Sudden singularities 

"Sudden singularities" are a recent type of cosmological milestones that have appeared 
in the literature. There could be past or future sudden singularities for which some 
time derivative of the scale factor diverges at finite time, while the scale factor itself 
remains finite. Again, future sudden singularities are more popular than past sudden 
singularities as those ones would be a most unusual and disturbing beginning to the 
history of the universe. 

Recently it has been speculated that in an expanding FLRW universe a curvature 
singularity may occur at a finite time before a "Big Crunch" for matter contents that 
satisfy both weak and strong energy conditions. This family of models has been fur- 
ther enlarged, and the same sort of behavior has also been found in inhomogeneous 
models Ej ■ 

We shall define the scale factor near a sudden singularity using the generalized 
power series expansion. 

Definition 4. Let the time of the sudden singularity, if one occurs, be t Q (past or 
future). A suitable definition of the exponent of a sudden singularity is to take r] = 
and rji > to give 

a(t) = c + d\t - t Q \ m + ... (3.48) 

with Co > and r/i non-integer. Thus a(t ) = Co is finite and a sufficient number of 
differentiations yields 

a {n) (t -> t Q ) ~ c ^(771 - l)(77i -2)...( Vl -n + l)\t- t Q | m ~" -> 00. (3.49) 
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The toy model considered by Barrow |Hl IE] can be written 

<t) = c [(t Q - tf - 1] + So(t - Uf (3.50) 

and falls into this classification when expanded around the time of the sudden singu- 
larity, t e , while it falls into the classification of big bang singularities considered above 
when expanded around the time of the big bang, t*. 

Extremality events 

Let's consider some other common cosmological milestones (which are not singularities 
in any sense) that are referred to as "extremality events". In particular, for these 
events, the scale factor a(t) exhibits a local extremum at finite time, and whence, 
a — > 0. In the vicinity of extremality events we can model the scale factor using 
ordinary Taylor series so that in terms of our generalized series we have r/o = and 

Vi e Z+. 

Definition 5. A "bounce" is any local minimum ofa(t), the time of such an event 
being denoted by t,, so that a^\t.) = 0. The "order" of the bounce is the first nonzero 
integer n for which the 2n 'th time derivative is strictly positive: 

a^ 2n \t.)>0, (3.51) 

so that 

a(t) = a(t.) + -^-a^ n \t.) [t - t.} 2n + .... (3.52) 

Note that a^(t.) = at the bounce is in agreement with H d 123 EH1 EH EDI 
EH! E21 EH E21 and that the 2n'th time derivative strictly positive is also mentioned 
in E21E|. 

Definition 6. A "turnaround" is any local maximum of a(t), the time of such an 
event being denoted by t , so that a^(t a ) = 0. The "order" of the turnaround is the 
first nonzero integer n for which the 2n 'th time derivative is strictly negative: 

a (2n) (t ) < 0, (3.53) 
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so that 

a(t) = a(t ) + -l-a^(t ) [t - t Q ] 2n + .... (3.54) 

Note that turnarounds are mentioned in El ESI • 

Definition 7. An "inflexion event" is an extremality event that is neither a local max- 
imum or a local minimum, the time of such an event being denoted by t@. (Inflexion 
events can be thought of as an extreme case of "loitering" , in the limit where the Hub- 
ble parameter momentarily vanishes at the inflexion event.) The order of the inflexion 
event is the first nonzero n for which 

a {2n+1 \t @ )^0, (3.55) 

so that 

a{t) = a{t @ ) + 7^TT7« (2n+1) (^) [t ~ t @ } 2n+1 + .... (3.56) 

The loitering universe scenario is an expanding Friedmann cosmology that under- 
goes a fairly recent phase of slow expansion. It is during this semi-static phase that 
large-scale structure would be formed. Loitering is also characterized by the fact that 
the Hubble parameter dips in value over a narrow redshift range referred to as the 
"loitering epoch" . During loitering, density perturbations would be expected to grow 
rapidly, and, since the expansion of the universe would slow down, its age near loitering 
would dramatically increase [HH1 EZI • 

Definition 8. The "emergent universe" of 131 \ can be thought of as an extremality 
event that has been pushed back into the infinite past. 

The "emergent universe" scenario consists of an inflationary universe that would 
emerge from a small static state that has within it the seeds for the development of 
the macroscopic universe. The universe would have a finite initial size, with a finite 
amount of inflation occurring over an infinite time in the past, and with inflation would 
then coming to an end via reheating in the standard way. The scale-factor would be 
bounded away from zero in the past and there would be no horizon problem and no 
singularity, since the initial state would be Einstein static. Also, the initial static state 
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could be chosen to have a radius above the Planck scale, so that these models could 
even avoid a quantum gravity regime |3*H I3U] . 

These definitions have been chosen to match with and simplify the definitions in 
articles jUj and [02] • Note that for bounces these definitions imply that there will be 
some open interval such that 

Vt G (t. - A, t.) U (t„ t. + A); a(t) > 0, (3.57) 

while for turnarounds there will be some open interval such that 

Vt G (t - A, t ) U (t D , t + A); a(t) < 0. (3.58) 

Note that unless the bounce or turnaround is of order one we cannot guarantee that 
at the extremality event itself a(t.) > or a(t D ) < 0. For inflexion events we can only 
assert the weaker condition of the existence of some open interval such that 

Vt G (t© — A, t@) U (t®, t© + A); d(t) has fixed sign. (3.59) 

Summary 

We have defined singular cosmological milestones (big bang, big crunch, big rip, sud- 
den singularity) and nonsingular cosmological milestones (extremality events) using a 
single framework based on generalized power series. From a kinematical point of view, 
major features will be deduced with the first two indicial exponents (7/ and 771). 

The value of the first indicial exponent 7/0 relates to the classification of the types 
of singular cosmological milestones as follows: 

• rjo > for big bangs or big crunches. 

• 7/o = for sudden singularities. 

• 7/0 < for big rips. 

For the nonsingular cosmological milestones (extremality events) , the parameteri- 
zation is simpler: only one positive integer is needed, i.e. the order of the extremality 
event, to determine the qualitative behaviour of the Taylor series of the scale factor 
a(t). 
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3.2.2 Spacetime curvature 

In this section, we will use the parameters define above to explore generic properties 
of the cosmological milestones. Do all cosmological milestones with scale factor singu- 
larities lead to curvature singularities? 



Cosmological parameters 

Before analyzing the spacetime curvature, we consider two cosmological parameters: 



a 



• the Hubble parameter H = — , characterizing the rate of expansion . 

a 

iii- a & tt-9 aa 

• the deceleration parameter 4 q = — H — -. 

a or 

Theorem 3. Consider a generalized power series expansion of the scale factor as 
defined in Definition QJ for rj ^ 0, that is for bangs, crunches, and rips, the Hubble 
parameter exhibits a generic l/(t — t & ) blow up, where t & is the the time of the event. 
Explicitly, 

lim H = / + °° 710 > ° ; (3.60) 

t^t+ — OO 7]o < 0. 

Proof. From the definition of a generic cosmological milestone 

a(t) = co\t - t & \ Vo + d\t - t & \ m + ... (3.61) 

we have (assuming t > t® for simplicity, otherwise one need merely formally reverse 
the flow of time) 5 

d(t) = c oV o(t - to)" " 1 + c m (t - to)" 1 " 1 + . . . (3.62) 



3 See section ETT1 
4 See section UTTl 

5 In fact for explicit calculations in the vicinity of any cosmological milestone it is always possible 
to choose the direction of time to force t > t© and so dispense with the need to take the absolute 
value \t — £q | , at least for one-sided calculations. This is not a physical restriction on the cosmological 
milestone, just a mathematical convenience which we shall adopt henceforth without further explicit 
discussion. 
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Keeping only the most dominant term, we have for rjo / 

a c r] (t - 1^ - 1 rjp 
a c {t - to^o t-t Q 
That is, for bangs, crunches, and rips the Hubble parameter exhibits a generic 
l/(t — t Q ) blow up. For T] > 0, lim t ^ tQ r] /(t — t Q ) = +oo whereas, for r] < 0, 
lim t ^ tQ 7] / (t - t & ) = -oo. □ 

ThcorcmElimplics that for the particular cosmological milestones such as big bangs, 
big crunches, and big rips, the Hubble parameter goes to infinity in the vicinity of the 
time of the event. 

Theorem 4. Consider a generalized power series expansion of the scale factor as 
defined in Definition^ for tjq = 0, that is either a sudden singularity or an extremality 
event, the Hubble parameter does not necessarily diverge at the cosmological milestone. 
As a matter of fact. 



7?i = 1; (3.64) 
m e (0,1); 



7] = 0; 

Jim H = <( ci/co Tj = 0; 

sign(ci)oo Tj = 0; 
where t Q is the the time of the event. 

Proof. For rj = 0, we need to go to the next highest term in the numerator (a term 
which depends on rji, which is guaranteed to be greater than zero by our definitions) 
to obtain 

h ~ ClT/l(t ~ t0)? ' 1 ' 1 =m-{t- to)*- 1 ; fob = 0;m>o). (3.65) 

Co Co 

In particular, this guarantees that power law behaviour is completely generic near 
the cosmological milestone, and while the value of the exponent is typically —1, there 
is an exceptional class of milestones (the sudden singularities and extremality events) 
for which the exponent will differ. It is not automatic that the Hubble parameter 
diverge at the cosmological milestone. In fact, 

i i?H = m £i (i _ ter - 1 = (° » = °; *>|; p.66) 

*-*s c [ a/ Cq 770 = 0; 771 = 1; 

□ 
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Theorem 0] implies that for the particular cosmological milestones such as sudden 
singularities and extremality events, the Hubble parameter does not necessary diverge 
and in particular, the Hubble parameter has a finite limit iff rj = 0, rji > 1, corre- 
sponding to a particular subset of the sudden singularities. 

The following theorem is a summary of the results encountered in Theorem |S1 and 
Theorem 0] 

Theorem 5. Summary: 

Consider a generalized power series expansion of the scale factor as defined in Defini- 
tion^ the Hubble parameter does not necessarily diverge at every different cosmological 
milestone. Indeed, 



lim H 



+oo 


ci/co 

sign(cx)oo 
— oo 



Vo > 

Vo = 

7/0 = 

r/ = 

Vo < 0, 



m > 1; 

Vi = 1; 
m e (o,i; 



(3.67) 



where t® is the the time of the event. 



Proof. Theorem El and Theorem 0] implies Theorem 



□ 



Now that we have analyzed the Hubble parameter in detail, we will consider the 

so-called deceleration parameter q = — H~ 2 = and give some similar results. 

a a 2 

Theorem 6. Consider a generalized power series expansion of the scale factor as 
defined in Definition QJ the deceleration parameter q can be either finite or infinite as 
the cosmological milestone is approached: 

• For bangs, crunches, and rips, (rjo ^ 0) the limit of q is always finite. 



For a certain subset of the sudden singularities (rj = 0) the limit q is infinite. 
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Indeed, 



(1 ~Vo)/vo 








sign(ci[l - r/i])oo 


7/0 = 


7/17^1; 







7/o = 


7/i = 1; 


7/2 > 2; 


-2c 2 cq/c\ 


r/o = 


T/i = 1; 


7/2 = 2; 


-sign(c 2 )oo 


7/0 = 


7/i = 1; 


7/2 e (1,2), 



lim q = { 7/ = 0; r/x = 1; r/ 2 > 2; (3.68) 



where t & is the the time of the event. 

Proof. First, we will consider the cosmological acceleration a: 

a{t) = c r/ (r/ - l)\t - t Q \^ 2 + — — t Q \^ 2 + c 2 r/ 2 (r/ 2 - l)\t - t Q \^~ 2 + ... 

(3.69) 

Remember that for the scale factor a(t), we have: 



a(t) = colt-tol" + c 1 |t-t | ?n +c 2 \t-t & \ m + c 3 \t-t Q \ m + 
Then provided rj ^ and 7/ 7^ 1 

a ^ 7/0(770 - 1) 
a ~ (t - £ ) 2 ' 

If ?7o = 0, then provided 771 ^ 1 

a 771(771 - l)ci 



c 



(t-t r~ 2 . 



If both 770 = and 771 = 1 



a 7/2(7/2 ~ l)c 2 , „ 2 _ 2 

- ~ (* - t Q ) • 

a c 



Finally, if 770 = 1, then (since 771 > 1) 

a 771(771 - l)ci 



Co 



(t-t )^- 3 . 



(3.70) 



(3.71; 



(3.72) 



(3.73) 



(3.74) 



The behaviour oi a/a near the milestone will be some power law, though the precise 
exponent of that power law will depend on the interplay between the various indicial 
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exponents fy. Note that there is at least one situation in which we have to calculate up 
to the third exponent rj 2 . We can now consider the so-called deceleration parameter 



q = -^ H - 2 = -^ (3-75) 
a a z 

Following the same sort of analysis, for r) ^ and r] ^ 1 we have the "generic" result: 

q~^—^; fa 7^0,1). (3.76) 
Vo 

For the "exceptional" cases we easily see: 

— If r/ = 0, then provided rji ^ 1 

fa -!)<*> 
77iCi 

— If both ?7o = and 771 = 1 

?~-^ife(i&-l)(f-f0r- 2 . (3-78) 
c i 

— Finally, if 770 = 1, then (since 771 > 1) 

m ( Vl - l) Cl {t _ t )m - K (3 . 79) 
c 

Again we see the ubiquity of power law behaviour, with a "generic" case and a limited 
number of special cases. Using the generic case and the special cases, we can calculate 
the limit of the deceleration parameter at the cosmological milestone. For ^ ^ and 
rjo ^ 1 we get the "generic" result: 

lim q = (1 - T) )/vo fa + 0, 1). (3.80) 

t— ><0 

For the "exceptional" cases we get: 

— If 7/0 = 0, then provided 771 ^ 1 

lim g = sign(ci[l — 771]) 00. (3.81) 
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If both r/o = and rji — 1 
lim q = 

Finally, if r/o = 1, then (since r/i > 1 



2c 2 c /c^ 
sign(c 2 )oo 



7/2 = 2; 
%G(1,2). 



lim o = 0. 

Note that in this case, the next dominant term r/i goes to as t 
therefore write this limit as 



lim q 



Vo)/Vo 



»70=1 



0. 



Hence, we have for the deceleration parameter: 



lim q 



(1 - r/o)/r/ 






sign(ci[l - 771] )oo 


Vo = 


117^1; 





r/ = 


m = i; 


-2c 2 cq/c\ 


Vo = 


m = i; 


-sign(c 2 )oo 


r/ = 


^1 = i; 



(3.82) 



(3.83) 



t Q . We can 



(3.84) 



V2 > 2; 
r/ 2 e(l,2) 



(3.85) 



Therefore, the deceleration parameter can be either finite or infinite, in particular, 
it has an infinite limit only for a certain subset of the sudden singularities r/ = 0, and 
for bangs and crunches the limit is always finite. This is largely because the definition 
of the deceleration parameter was carefully chosen to eliminate the leading t behaviour 
whenever possible. □ 



Polynomial curvature singularities 

Now to decide if the cosmological milestones we have defined are curvature singular- 
ities, we need to have a look at the Riemann tensor. The curvature is measured by 
the Riemann tensor from which can be constructed various scalar quantities, such as 
the Ricci scalar R, or higher order scalars R ab R a b, R abcd R a bcd, etc... If any of those 
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scalars go to infinity when approaching some point x, there is a polynomial curvature 
singularity at this point x. 

Because of the symmetries of FRW geometry, the Weyl tensor is automatically zero 
and so it suffices to consider the Ricci tensor (and implicitly the Ricci scalar). But 
because of spherical symmetry, and the perhaps less obvious translational symmetry, 
the only two non-zero orthonormal independent components of the Ricci tensor are 

R tt and R f f = R§§ = Rfo. (3.86) 

Thus to test for all possible polynomial curvature singularities it suffices to test for 
singularities in, for instance, R% and R?f. Alternatively one could consider R% and 
the Ricci scalar R, or even R^ and G^. 

From the metric I3~T| we can calculate the R% and Gn components in an orthonormal 
basis. We find: 

Rtt = -3 -; (3.87) 




(3.88) 



The interest of using those particular two combinations of the orthonormal com- 
ponents of the Ricci tensor is that they satisfy interesting properties: 

• they are linearly independent; 

• Rff is independent of the curvature of space (no k dependence); 

• Gn i s independent of a (minimizing the number of derivatives involved); 

• testing these two objects for finiteness is sufficient to completely characterize all 
polynomial curvature singularities in a FRW geometry. 

Theorem 7. The R% theorem: 

Consider a generalized power series expansion of the scale factor as defined in 
Definition QJ the R% component of the Ricci tensor in an orthonormal basis remains 
finite, as the cosmological milestone is approached, only provided: 
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• Vo = 0, 771 > 2; 



• Vo = 0, 771 = 1, 772 > 2; 

• 770 = 1 and rji > 3. 

T/iai is, except for these above rather limited cases, will blow up to infinity as the 
cosmological milestone is approached. 

Proof. The R% component of the Ricci tensor can be written as 

Rii = -3 ~ a (3.89) 

Using the generalized power series expansion we can rewrite the term only as a 
function of some (t — t Q ) factors in the vicinity of the event. 
Then provided r] 7^ and 770 7^ 1 

Rii =- 3 ~a~- 3 -(t-tvf- (3 - 90) 
If 7/0 = 0, then provided rji 7^ 1 



% = -3? ~ -3 ^ lfa " 1)Cl (t - t )—. (3.91] 
a c 



If both 770 = and 7/1 = 1 



ifc = -3^ ~ -3 m{V2 - 1)C2 (t - t )^ 2 . (3.92) 
a c 



Finally, if 770 = 1, then (since rji > 1) 



Rii = ~ -3 ^ lfa = 1)Cl (t - t )-- 3 . (3.93) 
a c 



The generic result is 



lim % = sign(r/ [l - r/ ]) 00 for 770 7^ 0; r/ 7^ 1, (3.94) 
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For the "exceptional" cases we get: 
- If r/o = 0, then provided 7/1 ^ 1 



lim i?.- 






r/o = 


= 0; 


m 


> 


2; 


-6ci/co 


r/o = 


= 0; 


r/i 




2; 


sign(ci)oo 


r/o = 


= 0; 


m 


e 


(1,2) 


— sign (ci) 00 


r/o = 


= 0; 


m 


e 


(0,1) 



If both 7/o = and 7/1 



lim 



Finally, if 7/0 = 1, then (since 7/1 > 1) 



lim i?.- 

t-t 1 



To summarize, in the vicinity of the cosmological milestone, we can write 



lim Rf 



—00 


r/o 


> 1 









r/o 


= 1 


r/i > 3; 




-18ci/c 


r/o 


= 1 


r/i = 3; 




— sign (ci) 00 


r/o 


= 1 


r/i G (0,3); 




+00 


r/o 


e (0,1); 







r/o 


= 


Vi >2; 




-6ci/c 


r/o 


= 


r/i = 2; 




sign (ci) 00 


r/o 


= 


r/i e (1,2); 







r/o 


= 


r/i = 1; 


r/2 


-6c 2 /c 


r/o 


= 


r/i = 1; 


r/2 


-sign(c 2 )oo 


r/o 


= 


r/i = 1; 


r/2 


— sign (ci) 00 


r/o 


= 


r/ie(0,l); 




—00 


r/o 


< 







(3.95) 






r/o 


= 0; 


'/1 


= 1; 


r/2 > 2; 




-6c 2 /c 


r/o 


= 0; 


r/i 


= 1; 


r/2 = 2; 


(3.96) 


-sign(c 2 )oo 


r/o 


= 0; 


r/i 


= 1; 


r/2 e (1,2); 








r/o = 


1; 


r/i > 3; 




-18ci/co 


r/o = 


1; 


r/i = 3; 


(3.97) 


— sign (ci) 00 


r/o = 


1; 


r/i G (0,3). 





(3.98) 



Therefore i?££ is finite provided 
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• Vo = 0, rji > 2; 



□ 



• Vo = 0,Vi = 1, 772 > 2; 

• 770 = 1 and 771 > 3. 

Theorem 8. The G% theorem: 

Consider a generalized power series expansion of the scale factor as defined in 
Definition QJ the G% component of the Ricci tensor in an orthonormal basis remains 
finite, as the cosmological milestone is approached, only provided: 



T]0 = ?7i > 1 

r/o = l 771 > 3 and Cq + k = 0. 



(3.99) 
(3.100) 



(Note that c\ + k = implies k = — 1, and Cq = \.) 

That is, except for these above rather limited cases, G% will blow up to infinity as 
the cosmological milestone is approached. 

Proof. The G% term can be written as a function of the scale factor, its first derivative 
and the curvature parameter k, 



a 2 k 
Git = 3 — H - 



\H 2 + f 2 



(3.10i; 



To analyze Gq recall that the condition that the Hubble parameter H remain finite 
was 770 = 0, 771 > 1 from Theorem In the vicinity of the milestone, we have 



G {{ ~ 3 



Vo 



k 



l(t-t Q 



c (t — £ 0< 



2% 



while 



Git ~ 3 



rfAit - t Q ) 2 ^~V + k 



Vo = 0. 



(3.102) 



(3.103) 



Whenever 770 > 1, the term k/c^{t — t & ) 2r, ° with k 7^ is dominant and it tends to 
00. In the case of 770 > 1 and k = 0, the term rfc/ (t — t & ) 2 is dominant and it tends 
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to oo. Whenever t] < 1 and 7] ^ 0, the term rfc/it — 1 ) 2 is dominant and therefore, 
G% 00 • For the special case where rjo = 1, is a difference of squares and one 



term can be balanced against the other. Finally, when 770 
be developed further using the indicial exponent rji. 
Therefore, after further calculations, we find, 



the first term needs to 



lim G t 



t~>tr. 



sign(fc)oo 


Vo 


> 1 








+00 


m 


> 1 


k = 0; 






sign(cQ + fc)oo 


m 


= 1 


cl + k ^ 0; 









m 


= 1 


c 2 + k = 0; 




>3; 


18C! 


m 


= 1 


c 2 + k = 0; 




= 3; 


sign(ci)oo 


m 


= 1 


c 2 + k = 0; 




G(l,3) 


+00 


Vo 


e (0,1); 






3k/c 2 


Vo 


= 


771 > i; 






3(c? + k)/cl 


Vo 


= 


Vi = l- 






+00 


Vo 


= 


m e (0,1). 







(3.104) 



Thus the necessary and sufficient conditions for G% to remain finite at the cosmological 
milestone are: 



Vo = 0, Vi > i; 

?7o = 1, cl + fc : 
(Note that cl - 



- 0, and 771 > 3. 
k = implies 



■1, and Co = 1.) 



□ 



Theorem 9. Polynomial curvature singularity theorem: 

Consider a generalized power series expansion of the scale factor as defined in 
Definition^ every cosmological milestone is a polynomial curvature singularity at the 
event except for the rather limited classes of cosmological milestones that satisfy the 
necessary and sufficient conditions: 

• T] = rji > 2 or; 

• 7] = rji = 1 T] 2 > 2 or; 
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• 7/o = 1 r]i > 3 and fc = — 1, c = 1. 

T/iai zs, /or a particular sub-class of sudden singularities (tjq = 0) and for a rather ex- 
ceptional type of big bang/crunch (as detailed above), there is no polynomial curvature 
singularity. 

Proof. From Theorem [7| and |H1 we have determined the conditions for Rn and Gq 
(respectively) to remain finite at the time the cosmological milestone occurs. 
For Rj£, those conditions are 

• ?7o = 0, % > 2; 

• Vo = 0, Vi = 1, V2 > 2; 

• 7/o = l and r/i > 3. 

For Gft, those conditions are 

• Vo = 0, 7/i > 1; 

• rjo — 1, c\ + A; = 0, and r/i > 3. 

(Note that Cg + fc = implies k = —1, and Co = 1.) 

Combining those two results, the necessary and sufficient conditions for both and 
Gn t° remain finite, so that a cosmological milestone is not a polynomial curvature 
singularity, are: 

• 7/o = 0, 7/1 > 2; 

• m = 0, 7/1 = 1, 7/2 > 2; 

• 7/0 = 1, = —1, c = 1, and 771 > 3. 

Note that the case 7/0 = 1 that corresponds to an exceptional type of big bang/ crunch is 
asymptotic to a Milne universe. Indeed, the special case a(t) = t, k = — 1 is called the 
Milne universe and is actually a disguised portion of Minkowski space. It corresponds 
to the interior of the future light cone based at some randomly specified point in 
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Minkowski space, with a spatial foliation denned by the proper time hyperboloids 
based on that event. Since this spacetime is flat, it corresponds to a universe which 
on the largest scales is empty. This is not a popular cosmological model. 

Any other conditions that are not mentioned above imply that and/or blow 
up at the time of the event; therefore, all other cases imply a polynomial curvature 
singularity. □ 

Even if there are a few classes of cosmological milestones that are not polyno- 
mial singularities, one can wonder what happens to them when looking at finite-order 
derivatives of the curvature tensor. Is it still possible to find an interesting type of 
singularity where the derivative curvature singularity does not blow up? 

Derivative curvature singularities 

Definition 9. A derivative curvature singularity is defined by some polynomial con- 
structed from finite-order derivatives of the curvature tensor blowing up. 

In our case, because of the symmetries of the FRW universe the only interesting 
derivatives will be time derivatives, and so the only objects we need to consider are 

^ and (3.105) 
d n t d n t v ' 

Theorem 10. Derivative curvature singularity theorem 

The only two situations in which a cosmological milestone is not a derivative cur- 
vature singularity are if: 

• 7] = 0, rji G Z + ; corresponding to an extremality event (bounce, turnaround, or 
inflexion event) rather than a bang, crunch, rip, or sudden singularity; 

• i]q = 1, k = —1, c = 1, r]i G Z + , and rji > 3; corresponding to a FRW geometry 
that smoothly asymptotes near the cosmological milestone to the Riemann-flat 
Milne universe. 

Proof. We can write d n R^/d n t and d n G^/d n t using the scale factor a(t) and its deriva- 
tives oW(t), 

d n Rn a(' ra+2 ) , , , . 

= —3 h (lower-order derivatives); (3.106) 

d n t a 
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while 

d n G~ a a (n+1) 

— - — — = 3 h (lower-order derivatives). (3.107) 

d n t a 2 

Thus to avoid a rath-order derivative curvature singularity we must at the very 
least keep a ( - n+2 - ) / a finite, and furthermore all related lower-order derivatives of the 
form a^/a, with j < n + 2, must also be finite. To prevent any arbitrary-order 
derivative singularity from occurring, that is for both d n R^/d n t and d n G^/d n t to 
remain finite for all n G Z + , we must force all a^/a to remain finite. This condition 
holds in addition to the constraint coming from polynomial curvature singularities 
derived above. 

The first constraint to hold, so that there is no polynomial curvature singularity, 
is rjo = 0. In this case, the scale factor can be written as: 

a(t) = c + Cl \t - + c 2 \t - t & \ m + c 3 \t - t e \ m + ... (3.108) 

Remember that the dominant term in the scale factor is Cq which implies that a(t) 
finite at the time of the event. Therefore, for /a to remain finite it suffices that 
be finite. 

The second constraint to hold, so that there is no polynomial curvature singularity, 
is rjo = 1. In this case, the scale factor can be written as: 

a(t) = c |t - * | + c x \t - to)" 1 + c 2 \t - t & \^ 2 + c 3 \t - to]" 3 + . . . (3.109) 

Remember that the dominant term in the scale factor is Co\t — t & \ which implies that 
a(t) is not finite at the time of the event. However, the term \t — from a(t) can be 
"absorbed" in the term a^ n+2 ^ as follows: 



a 
a 







E, 


dTji . . 


■ (Vi - J + 1) 1 


t-t Q \> 


n-i 




t-t Q 


+ Ci 


t-t Q 


171 + c 2 


t-t Q 


\n* + c 3 \t-t Q 


m _|_ , . , 



(3.110) 



keeping only dominant terms, we have 



Ei °iVi ■■■{Vi- 


-J + l 


) 


t-t e 


ni 3 E l <kVi ■■■(Vi-j + 1 ) 


t-t Q 


W-3- 1 


a c 


t-t Q 


Co 



(3.111; 

Therefore, for /a to remain finite it suffices that be finite. 
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These two constraints also imply that to force all a^' to remain finite, we must force 
all the indicial exponents rji to be non-negative integers, thus making a(t) a Taylor 
series. Indeed, we have 



a 



= CiTji . . . {rji - j + 1) \t - t, 



■0 



(3.112) 



and is finite if and only if 




finite, 



(3.113) 




finite: otherwise, there will be eventually a value of j that will make the term rji—j < 0, 
to avoid that rji must be integers so the term \t — t & \ m ~^ disappear and do not tend 



Hence, we have demonstrated that almost all cosmo logical milestones are physical 
singularities, apart from a very limited sub-class corresponding to either extremality 
events (bounces, turnarounds, inflexion events) or an asymptotically empty universe. 

3.3 Cosmological milestones and dynamics. 

In the previous section, we have considered only kinematics, from now on, we will 
use the Friedmann equations and the Einstein equations to take our definitions to a 
dynamic level. 

3.3.1 Introduction 

To now start to include dynamics we relate the geometry to the density and pressure 
using the Friedmann equations and then ask what happens to the energy conditions 
at the cosmological milestones. Even though we may have good reason to suspect 
that the energy conditions are not truly fundamental [3] : they make a very good first 



to infinity. 



□ 
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pass at the problem of quantifying just how "strange" physics gets at the cosmological 
milestone. 



Using the Friedmann equations, we can write the pressure p(t) and the density p(t) 
using the scale factor and its first two derivatives: 

= 3 (M)' < 3 - 115 > 

, , Ob CL Jv , s 

pit) = -2 3.116 

p(t) + 3p(t) = -6-. (3.117) 

a 



Therefore, we have 



and 



a a z 



(d a 2 k \ , 
- + 2- + 2- . (3.119) 
a a z a z J 

The standard energy conditions are the null, weak, strong, and dominant energy 
conditions 6 which for a FRW spacetime specialise to [13 0] : 

• [NEC] p + p > 0. 

In view of equation (|3. 118(1 this reduces to 

a < ; that is k > a a — a . (3.120) 

a 

• [WEC] This specializes to the NEC plus p > 0. 
This reduces to the NEC plus the condition 



a 1 + k > 0; that is k > -a 2 . (3.12i; 



[SEC] This specializes to the NEC plus p + 3p > 0. 
This reduces to the NEC plus the deceleration condition 



a < 0. (3.122) 



3 See section T2.6I for more details on the energy conditions. 



57 



• [DEC] p±p > 0. 

This reduces to the NEC plus the condition 

2(d 2 + k) , , (a a + 2d 2 ) 
d> ^— -; that is fc > — - -. (3.123) 

The null energy condition NEC is the most interesting condition because it is the 
weakest of the standard energy conditions and it leads to the strongest theorems. 

3.3.2 NEC 

We have seen that the NEC is satisfied iff k ^ a a — a 2 . 

Theorem 11. The NEC is definitely satisfied (meaning the inequality is strict ) at a 
generic cosmological milestone iff: 

• Big bangs/crunches: 

— T] > I, k = +1; 

— 770 = 1: for k — +1 and k = 0, and also for k = — 1 with the proviso that 
co > 1 ; 

— 7/0 G (0, 1): for any value of k; 

• Sudden singularities /extremality events: 

~ Vo — subject to the additional constraints: 

* r/i > 2 and k = +1; 

* 771 = 2 and fc > 2c Ci/ 

* r)i e (1, 2) and Ci < 0; 

* 771 = 1, 7/2 > 2: /or fc = +1, k = 0, and /or k — — 1 win ide proviso 
\ci\ > 1; 

* r/i = 1, 7/2 = 2: /or fc > 2c c 2 - c 2 ; 

* 7/1 = 1, 7/2 G (1, 2), and c 2 < 0; 

* rj 1 G (0, 1), and Ci > 0. 
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Theorem 12. The NEC is marginally satisfied (meaning the non-strict inequality is 
actually an equality) at a generic cosmological milestone iff: 

• Big bangs/crunches: 

— 770 > 1, k = 0; 

— rjo — 1, k = — 1 with the proviso that cq — 1; 

• Sudden singularities /extremality events: 

— rjo = subject to the additional constraints: 

* r}\ > 2 and k = 0; 

* rji = 2 and k = 2c^C\ (which requires k 7^ 0); 

* rji — 1, r]2 > 2: for k = — 1 and c\ = ±1; 

* ?7i = 1, 772 = 2, and k = 2c c 2 — c\. 

Theorem 13. The NEC is definitely violated (the inequality is strictly violated,) at 
a generic milestone iff: 

• Big bangs/crunches: 

— 770 > 1, k = -I; 

— 770 = 1 : /or = — 1 u>it/j t/ie proviso that cq < 1 ; 

• Sudden singularities /extremality events: 

— rjo = subject to the additional constraints: 

* rji > 2 and k = —1; 

* rji = 2 and < 2coCi/ 

* 771 G (1, 2) and c\ < 0; 

* rji = 1, rj2 > 2: /or fc = — 1 u>it/i t/ie proviso \ci\ < 1; 

* rj\ = \,rj 2 = 2: for k < 2c c 2 - cf; 

* 771 = 1, 772 6 (1, 2), and c 2 > 0; 
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* rji E (0, 1), and c\ < 0. 
• Big rips: 
- r/ < 0. 

In particular, "big rips" (t/q < 0) will violate the NEC in the vicinity of the milestone 
and therefore, they will also violate every other energy conditions. 

Proof. The NEC is satisfied if k ^ a a — a 2 . Remember that near the cosmological 
milestone we can write the scale factor as a generalized power series as defined in 
Defintion ^ 

a(t) = co|t - t r° + cx\t - t Q \ ni + c 2 \t - t Q \ 

Therefore, near any generic cosmological milestone 

2 



|i/i> 



+ C 3 \t - t Q \ r " + 



a a — ' 
while in the degenerate cases 



and 



a a - a 2 ~ c ci 771(771 - 1) (t - t Q ) {m 2) ; 



a a - a 2 ~ c c 2 772(772 — 1) (t — t©) (r?2 2) 



(Vo + 0), 
(t/o = 0;?7i i- 1), 
(Vo = 0; 771 = 



(3.124) 
(3.125) 
(3.126) 
(3.127) 



After further calculations, we find the limit of (a a — a 2 ) when the milestone occurs: 



lim (a a 
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—00 


?/o e (0,1) 
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2Q)Ci 
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+sign(ci)oo 
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2CoC2 — c\ 
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+sign(c 2 )oo 
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7/1 = 1; 
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— sign (ci) 00 
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7/1 e (0,1) 




+00 
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(3.128) 
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The NEC is definitely satisfied iff: 

k>aa-a 2 , (3.129) 

that is, iff: 

• rjo > 1, k = +1; 

• ?7o = 1: for = +1 and k — 0, and also for k = —1 with the proviso that Co > 1; 

• 7/ G (0, 1): for any value of fc; 

• ?7o = subject to the additional constraints: 

— r]i > 2 and k = +1; 

— 7/1 = 2 and > 2cqCi, 

— 7/i G (1, 2) and c x < 0; 

— 7/i = 1, 7/2 > 2: for fc = +1, k — 0, and for fc = — 1 with the proviso |ci| > 1; 
-7/1 = 1,7/2 = 2: for fc > 2c c 2 - cj; 

— rfi = 1, rj 2 G (1, 2), and c 2 < 0; 

— r/i G (0, 1), and ci > 0. 

This proves Theorem II II 

The NEC is marginally satisfied iff: 

k = a a — d 2 , (3.130) 

that is, iff 

• T] > 1, k = 0; 

• 7/o = l, k = —1 with the proviso that Co = 1; 

• r/ = subject to the additional constraints: 

— 7/i > 2 and fc = 0; 
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— r/i = 2 and k = 2c Ci (which requires k ^ 0); 

— 7/i = 1, 7/2 > 2: for = — 1 and C\ = ±1; 
-7/1 = 1,7/2 = 2, and = 2c c 2 - cf. 

This proves Theorem H"2l 

The NEC is definitely violated iff: 

k<aa-a 2 , (3.131) 

that is, iff 

• 7/o > 1, k = -1; 

• 7/o = l: for = — 1 with the proviso that cq < 1; 

• 7/o = subject to the additional constraints: 

— 7/1 > 2 and k = —1; 

— 7/i = 2 and & < 2coCi; 

— 7/1 G (1, 2) and Ci < 0; 

— 7/1 = 1, 7/2 > 2: for = — 1 with the proviso |ci| < 1; 
-7/1 = 1,7/2 = 2: for fc < 2c c 2 - c\; 

— 7/1 = 1, 7/2 G (1, 2), and c 2 > 0; 

— 7/1 G (0, 1), and c\ < 0. 

• 7/o < 0. 

This proves Theorem ITBI □ 

Theorems ^2 El an d El present a pattern that is in agreement with the stan- 
dard folklore, and is a systematic expression of results that are otherwise scattered 
throughout the literature. For instance, it is immediate that "big rips" always vi- 
olate the NEC (and hence all the other energy conditions), and that for big bangs 
and crunches the range 7/0 G (0, 1) is preferred, since it is only in this range that the 
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NEC holds independent of the spatial curvature. In particular, if the NEC is to hold 
independent of the sign of space curvature k all the way down to the singularity, then 
in the vicinity of the singularity the dominant term in the scale factor a(t) is bounded 
from both above and below by 

Co(t - U) < adominant(t) < C . (3.132) 

Note that for sufficiently violent big bangs (770 > 1) and hyperbolic spatial curvature 
the NEC is violated: this indicates that "phantom matter" need not always lead to 
a "big rip" ^B]; it might in contrast lead to a particularly violent bang or crunch. It 
is for ?7o = (which corresponds to either sudden singularities or extremality events) 
that the analysis becomes somewhat tedious. Certainly, many (though not all) of the 
sudden singularities violate the NEC. 

3.3.3 WEC 

We have seen that the WEC is satisfied iff 

• the NEC is satisfied, that is if k ^ a a — a 2 , 

• and, p ^ 0, that is if k ^ —a 2 . 

Theorem 14. The WEC is definitely satisfied (meaning the inequality is strict ) at a 
generic cosmological milestone iff the NEC is definitely satisfied (Theorem \ll\) . 

Theorem 15. The WEC is marginally satisfied (meaning the non-strict inequality 
is actually an equality) at a generic cosmological milestone iff the NEC is marginally 
satisfied (Theorem X 1%) . 

Theorem 16. The WEC is definitely violated (meaning the inequality is strict vi- 
olated ) at a generic cosmological milestone for any cases that are not mentioned in 
Theorems^^ or\l 6 A 

Proof. The WEC is satisfied if k ^ a a — a 2 and k ^ —a 2 . The first condition 
k^a'd — a 2 has been dealt with in the proof of Theorems II 1| Upland 1131 We need to 
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look at k ^ — d 2 . Near the cosmo logical milestone we can write the scale factor as a 
generalized power series as defined in Definition^ 



a(t) = c \t - t^o 



I '/i 



c 2 \t - t 



0l 



c 3 \t-t 



Ci\t — t Q 

Therefore, near any generic cosmological milestone 



0l 



-d 2 ~ 



while in the degenerate case 



—a ~ 



-^(t-te) 2 ^" 1 ); (770 = 0). 



(3.133) 



(3.134) 



(3.135) 



Therefore 
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— oo 
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(3.136) 



Since we want the condition fc > — d 2 to hold in addition to the NEC, it follows that 
the this condition is equivalent to the NEC for bangs crunches and rips, and that 
the only changes arise when dealing with 770 = (sudden singularities or extremality 
events), but overall the conditions on the NEC are dominant. Therefore, the WEC 
holds if and only if the NEC holds. 

This proves Theorems El El and El D 

Certainly, many (though not all) of the sudden singularities violate the WEC. 
In particular, if the WEC is to hold independent of the sign of space curvature k all 
the way down to the singularity, then in the vicinity of the singularity the dominant 
term in the scale factor a(t) is bounded from both above and below by 



Colt — t*) < adominant(^) < Q)- 



(3.137) 
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3.3.4 SEC 

We have seen that the SEC is satisfied iff 

• the NEC is satisfied, that is if k ^ a a — a 2 , 

• and, p + 3p ^ 0, that is if a ^ 0. 

Theorem 17. The SEC is definitely satisfied (meaning the inequality is strict,) at a 
generic cosmological milestone iff: 

• Big bangs/crunches: 

— f]o = 1: for k = +1 and k = 0; 

— f]o £ (0, 1): for any value of k; 

• Sudden singularities/ extremality events: 

— ?7o = subject to the additional constraints: 

* rji > 2, k — +1 and c\ < 0; 

* rji = 2, k > 2c Ci and C\ < 0; 

* r/x G (1, 2) and c\ < 0; 

* r]i = 1, f]2 > 2 and C2 < 0: for k = +1, k = 0, and for k = — 1 with the 
proviso |ci| > 1; 

* i]i = 1, r/2 = 2 and Ci < 0: for k > 2c C2 — c\; 

* 771 = 1, 772 G (1, 2), and c 2 < 0; 

* ()i G (0, 1), and c\ > 0. 

Theorem 18. The SEC is marginally satisfied (meaning the non-strict inequality is 
actually an equality) at a generic cosmological milestone iff: 

• Big bangs/crunches: 

— 770 = 1; Co = 1, C\ < 0, and k = — 1; 

• Sudden singularities/ extremality events: 
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7/o = subject to the additional constraints: 

* 771 > 2, k — and ci < 0; 

* r]\ — 2, k — 2cqC\ and c\ < 0; 
(which requires that Co = —l/{2c\)) 

* 7/1 = 1, 7/ 2 > 2 and c 2 < 0, fc = —1, C\ = ±1 and c 2 < 0; 

* 7/1 = 1, 7/2 = 2 and c 2 < 0, and k = 2c C2 — cf; 

(which requires that k = — 1 and therefore 2cqC2 — c 2 = - 



Theorem 19. TTie SEC is definitely violated (meaning the inequality is strict vio- 
lated,) at a generic cosmological milestone for any cases that are not mentioned in 
Theorems^^ or\18l 

Proof. The SEC is satisfied if k ^ a a — a 2 and a ^ 0. The first condition k^a'd — a 2 
has been dealt with in the proof of Theorems ^2 EH an d ED We need to look at d ^ 0. 
Near the cosmological milestone we can write the scale factor as a generalized power 
series as defined in DefintionQ 



+ c 2 \t - t Q \" 2 + c 3 \t - t Q \" 3 + 



a(t) = co|* — t&\ V0 + ci|t - t e 
Therefore, near any generic cosmological milestone 

d~ 7/0(7/0-1) c t ( "°- 2) ; (t/o^O, T/o^l) 
while in the degenerate cases 



0; 
1: 



"17^1; 



772(772 - 1) c 2 (t 



Therefore 



lim [sign (a)] 



sign(ci) 
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sign (7/1 [t/1 
sign(c 2 ) 
+ 1 



(vo = 0; ^i = l;), 

% > 1 

Vo = l; 

% e (0, 1) : 

77o = 0; 7/1 ^ 1; 

?/o = 0; 7/1 = 1; 

Vo < 0. 



(3.138) 
(3.139) 

(3.140) 
(3.141) 



(3.142) 
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Since we want the condition a < to hold in addition to the NEC, it follows that the 
SEC is definitely satisfied when those conditions below hold: 

• ?7o = 1: for k = +1 and k = 0; 

• rjo G (0, 1): for any value of k; 

• rjo = subject to the additional constraints: 

— rji > 2, k = +1 and c\ < 0; 

— rji — 2, k > 2c ci and c\ < 0; 

— T)\ G (1, 2) and c x < 0; 

— 771 = 1, 772 > 2 and C2 < 0: for = +1, k — 0, and for /c = — 1 with the 
proviso |ci| > 1; 

— 771 = 1, 772 = 2 and c 2 < 0: for fc > 2c c 2 — c^; 

— 771 = 1, 772 G (1, 2), and c 2 < 0; 

— 771 G (0, 1), and c x > 0. 

The SEC is marginally satisfied when: 

• ?7o = 1, Co = 1, C\ < 0, and k = — 1; 

• ?7o = subject to the additional constraints: 

— 7/1 > 2, = and ci < 0; 

— r/i = 2, k = 2cqC\ and ci < 0; 
(which requires that c = — l/(2ci)) 

— 771 = 1, 772 > 2 and c 2 < 0, fc = —1, c\ = ±1 and c 2 < 0; 

— 771 = 1, 772 = 2 and c 2 < 0, and = 2c c 2 — cf; 

(which requires that k = — 1 and therefore 2coc 2 — c\ = —1) 

These conditions prove Theorems \H\ [HD and El D 
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The SEC is definitely violated for all "rips", and is also definitely violated for all 
"violent" bangs and crunches with rjo > 1. The SEC is definitely satisfied for bangs and 
crunches in the range t] G (0, 1). In particular, if the SEC is to hold independent of 
the sign of space curvature k all the way down to the singularity, then in the vicinity 
of the singularity the dominant term in the scale factor a(t) is bounded from both 
above and below by 

Co{t - U) < a d ominant(^) < C . (3.143) 

For turnarounds or bounces r\\ = 2n, n G Z + , and the sign of C\ governs possible 
SEC violations: Turnarounds satisfy the SEC while bounces violate the SEC 7 . For 
inflexion events, where r\\ — 2n + 1, n G Z + , the SEC is violated either just before or 
just after the inflexion event 8 . 



3.3.5 DEC 

We have seen that the DEC is satisfied iff 

• the NEC is satisfied, that is if k ^ a a — a 2 , 

(ad + 2d 2 ) 

• and, p — p ^ 0, that is if k > . 

Theorem 20. The DECis definitely satisfied (meaning the inequality is strict ) at a 
generic cosmological milestone iff: 

• Big bangs/crunches: 

- Tjo > 1, k = +1. 

— 770 = 1, k — 0, +1, and k = — 1 if Co > 1. 

— 770 G (1/3, 1), for all k. 

- Tjo = 1/3, 

* rji > 5/3, k = +1. 

7 See also 03 
8 See also [H|E2]. 
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* ?7i = 5/3, k > -fcocx. 

* rn E (1/3, 5/3), ci > 0, for all k. 

• Sudden singularities /extremality events: 

- Vo = 0, 

* 771 > 2, k = +1. 

* 771 = 2, > Co max{2ci, — ci}. 

* rji = 1, 772 > 2, k = 0, +1, and = — 1 |ci| > 1. 

* 771 = 1, 772 = 2 ; A; > Co max{2ci, — ci} — cf . 

Theorem 21. TTie £?£'C is marginally satisfied (meaning the non-strict inequality is 
actually an equality) at a generic cosmological milestone iff: 

• Big bangs/crunches: 

- r) > 1, k = 0. 

- rj = 1, k = -1 if c = I. 

- rjo = 1/3, 

* 771 > 5/3, k = 0. 

* 771 = 5/3, k = —cqCx. 

• Sudden singularities /extremality events: 

- Vo = 0, 

* 771 > 2, k = 0. 

* rji = 2, wift k = cq max{2ci, — Ci}. 

* m = 1, m > 2, k = -1 i/ |ci| = 1. 

* rji — 1 } r)2 — 2, k — cq max{2ci, — c{\ — c\. 

Theorem 22. The DEC is definitely violated (meaning the inequality is strict vi- 
olated ) at a generic cosmological milestone for any cases that are not mentioned in 
Theorems\EB orWR 



69 



Proof. The DEC is satisfied if k ^ a a — a 2 and k > —(a d + 2d 2 )/2. The first condition 
k^aa — a? has been dealt with in the proof of Theorems ^2 El an d ED We need to 
look at k > —(a a + 2d 2 )/2. Near the cosmological milestone we can write the scale 
factor as a generalized power series as defined in Defintion Q 

a(t) = co|t - t r° + cx|t - t Q \^ +c 2 \t- t |" 2 +c 3 \t-t Q \^ + ... (3.144) 

Therefore, near any generic cosmological milestone 

_(aa + 2a a ^_ , (3^-l)cg (t _ ts)2(w _ 1); ( ^ Q ^ ^ ^ ^ 
while in the various degenerate cases 

_(«a + 2^ ^ - 1) ^ c (f _ ie)( „-2, ho = . (X146) 

_ (aa + 2g) ^ _ »(,» -1) Co c 2 ft _ tB)fe _ 2) _ ^ ()fc = 0; |h = l)i (3147) 

_ (aa + 2a') ^ _ (3„, + 2)(3,,. - 1) c„ c, (f _ ^fr,,^. ta=1/3) . (3 . 148) 
2 18 
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Therefore 
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+oo 




'/O 


e (0,1/3); 











'/o 


= 
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= 




m = 


2; 




-sign(ci 


)oo 


'/o 


= 
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(3.149) 

Remember that to satisfy the DEC one needs to satisfy the NEC in addition to the 
constraint coming from the above. Let us define 



K = max < lim lad — a 2 ] , lim 



(ad + 2d 2 



(3.150) 



Then satisfying the DEC at the cosmological milestone is equivalent to the constraint 



k > K. 



(3.15i; 
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For bangs and crunches we calculate: 





K 



-oo 

14 

--gCod 

-sign(ci)oo 
-oo 



For rips we have: 



K = +00 



m > l; 
??o = l; 
770 e (1/3,1); 
7/0 = 1/3 
Vo = 1/3 
7/0 = 1/3 
%G (0,1/3). 



770 < 0. 



Vi > 5/3 
r/i = 5/3 
Vi < 5/3 



(3.152) 



(3.153) 



Note that rips, because they violate the NEC, always violate the DEC, so the particular 
constraint derived above is not the controlling feature. 
— Finally, for sudden singularities and extremality events: 



K 








Vo 


= 


Vi > 2; 




Co max{2ci, 


-ci} 


Vo 


= 


r/i = 2; 




+00 




Vo 


= 


me (1,2) 




-cl 




Vo 


= 


vi = i; 


V2 > 2; 


c max{2ci, 


-ci} - cl 


Vo 


= 


vi = i; 


r/ 2 = 2; 


+00 




Vo 


= 


vi = i; 


r/ 2 G(l,2 


+00 




Vo 


= 


771 e (0,1). 





(3.154) 



As we have now seen (several times), the r/ = case is "special" and requires extra care 
and delicacy in the analysis — this is the underlying reason why "sudden singularities" 
are so "fragile" , and so dependent on the specific details of the particular model. Indeed 
there are several classes of cosmological milestone for which the DEC is satisfied. For 
instance, a complete catalogue of the bangs and crunches for which the DEC is satisfied 
is: 

• r)o>l,k = 0,+l. 

• t]q = 1, k — 0, +1, and k = — 1 if c > 1. 
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• rio G (1/3, 1), for all k. 

• T] = 1/3, 77x > 5/3, = 0,+l. 

• Tj = 1/3, 771 = 5/3, k > —cad. 

• 770 = 1/3, 771 G (1/3,5/3), ci > 0, for all k. 

Similarly, a complete catalogue of the sudden singularities for which the DEC is 
satisfied is: 

• ?7o = 0, 77! > 2, k = 0,+l. 

• rjo — 0, rji — 2, with k > Cq max{2ci, — ci}. 

(In particular this requires k = +1 though that is not sufficient.) 

• % = 0, r)x = 1, 772 > 2, k = 0, +1, and k = — 1 if |ci| > 1. 

• 77o = 0, 771 = 1, 772 = 2, fc > c max{2ci, -cj - cf . 

These catalogues prove Theorem EDI EH and E3 □ 

If one wishes to use "normal" matter (that is, matter satisfying all of the energy 
conditions) to drive a bang or a crunch independent of the sign of space curvature k, 
then one is forced in a model independent manner into the range 770 G (1/3,1). In 
particular, if the DEC is to hold independent of the sign of space curvature k all the 
way down to the singularity, then in the vicinity of the singularity the dominant term 
in the scale factor a(t) is bounded from both above and below by 

Co(t - £*) < GtdominantW < C (t - t*) 1/3 . (3.155) 

The fact that we find special cases of sudden singularities that do not violate the 
DEC appears at first glance, to contradict the analysis of Lake [52] , who claimed that 
all sudden singularities violate the DEC. The resolution of this apparent contradiction 
lies in the question "just how sudden is the sudden singularity?" Lake takes as his 
definition d(t — > t & ) = — oo, corresponding in our analysis to 770 = with either 
771 G (0, 1) U (1, 2) or 771 = 1 with 7/2 < 2. In this situation our results certainly agree 
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with those of Lake: if a(t — > t & ) = — oo then the DEC is certainly violated. However, 
our "counterexamples" where the DEC is satisfied all satisfy d(t — > t Q ) = finite, and 
for these counterexamples it is only some higher derivative cS n \t — > t Q ) for n > 3 that 
diverges. Sudden singularities of this type are sufficiently "gentle" that at least for 
some of them even the DEC can be satisfied all the way to the singularity: a conclusion 
completely in agreement with Barrow and Tsagas [H], though now we have a complete 
characterization of those situations for which the DEC can be satisfied. 



3.4 Total age of the universe 

In this section, we express the density as a function of the scale factor p(a(t)) and 
analyze the Friedmann equation ()3.2j) rewritten with p(a(t)). This equation can then 
be formally solved by integration to find t(a), time as a function of the scale parameter. 

We adopt a barotropic equation of state by assuming the statement that the pressure 
is a function of density, p(p). If we express the density as a function of the scale factor, 
we can write p = p(a(t)). We can then write, 

T- = ( 3 - 156 ) 
da da I at 

Now rewriting the conservation equation, in units where 8tiGn = 1 and c = 1, we 
obtain, 

~p+p(pY 



da 



so that 



Integrating, we get 



Definition 10. Define, 



0, (3.157) 

dp da 



3[p + p(p)] a' 
dp f da 



3[p + p(p)] J a 



(3.158) 
(3.159) 



h(p) = j dp . (3.160) 

' [P + P(P)} 
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where p p= o is the density at zero pressure which is defined by 

p(p P =o) = 0. (3.161) 

The function h(p) is called the relativistic specific enthalpy. 

Using Definition EH then in terms of the present day density po and the present 
day scale factor a 

-h(p)-h( Po ) = \n(a/a ). (3.162) 



We can now write 

1 

a = Qq exp 



3 (%)-% )) 



(3.163) 



and deduce the existence of a function / such as 

a = f(p). (3.164) 

Note that when the density increases, assuming the NEC holds so does the enthalpy 
and therefore, the scale factor decreases. Now, apply the inverse function theorem to 
the function / in equation (j3.164j) . After all 

dh 1 , 

— = — ^ 3.165 

dp p + p{p) 

except for very peculiar isolated cases. Therefore, the conditions for employing the 
inverse function theorem apply, and thus we can write 

p = r\a) (3.166) 



which we often just abbreviate by writing p(a). 
The Friedmann equation can now be written as: 



p(°(*)) = 3('^ + -£d 



a 2 (t) a 2 (£). 

Equation (|3.167|) can be rearranged into equation (J3.168)) and formally solved by 
integration to find t(a), time as a function of scale parameter a. First note 

-p(a(t)) a 2 {t) - d 2 {t) = k. (3.168) 
3 
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Theorem 23. Let the time now be t and the size of the universe now be do, the 
standard solution of equation \3. 1 68\) in the current epoch is: 

t(a)= f da ^ . + tQ ( 3169 ) 

\ \p{a)d 2 — k 

Proof. There are three solutions to equation ()3.168|) : 
1. 

a(t) = p(a(t)) a 2 (t)/3 = k, (3.170) 

this means that the size of the universe does not change over time. This case 
is the Einstein static universe. A more complicated analysis would show that 
for "normal" matter the Einstein static universe is always unstable. It might 
be possible to find some form of "unusual" matter to make the Einstein static 



universe stable. 



2. 



in this case 



t(a) — — I — -. V constant, (3.171) 

J^p{d)a 2 — k 

f a da . 

> =>• / > 0. (3.172) 



|p(a)a 2 — k J J |p(a)a 2 — k 



Time decreases as the integral increases. That is time increases as the integral 
decreases (as the scale factor a decreases). So time increases as the universe gets 
smaller. This represents a contracting universe. Our universe is not currently 
contracting, but it is conceivable that the cosmological expansion could achieve 
a maximum, then turn around, and start to contract. If this happens then we 
would want to use this solution during the contracting phase. 

t(a) = [ — . = h constant, (3.173) 

J J |p(a)a 2 — k 

this solution is the standard solution appropriate to the current cosmological 
epoch. 



76 



Therefore, if we define t as the time now and ao as the size of the universe now, the 
solution of equation (|3.168|) is: 

t(a)= [ da _ + 1 (3.174) 

V |p( a ) fl2 — k 

□ 

Domains of reality 

We have shown that the current standard solution to equation ()3.168|) contains an 
integral of a square root \J p(a)a 2 /3 — k. However, this integral exists and is real only 
if the term inside the square root is positive. 



There are at least 5 situations wherein the square root y p(a)a 2 /3 — k is real: 

1. The square root is real only on [0,a max ]. (This corresponds to a big bang, and 
then a turnaround and a possible recollapse, or a future sudden singularity.) 

2. The square root is real only on [a*, oo]. (This corresponds to a bounce or a past 
sudden singularity then followed by infinite expansion.) 

3. The square root is real on the entire positve axis [0, oo]. (This is the most 
common assumption in cosmology, it corresponds to a big bang, then followed 
by infinite expansion). 

4. The square root is real only on a finite interval bounded away from zero, [a*, a max ]. 
(This corresponds to a bounce or past sudden singularity subsequently followed 
by a turnaround or future sudden singularity; a possible oscillating universe...) 

5. Disjoint unions of various of the above. (Not physically interesting.) 
3.4.1 Age of the universe r 

In this section, we will concentrate on the lower bound of the intervals of the different 
cases mentioned above, where the square root y 7 p(a)a 2 /3 — k is real. But first, we 
will present a generic definition and a generic theorem on the age of the universe. 
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Definition 11. The age of the universe r is the time elapsed since a = a . 

Here a Q can represent a big bang (in this case a© = 0) or a minimum value (for a 
bounce or past sudden singularity in this case a Q = a*). 



Theorem 24. The age of the universe r is given by, 

rao da 



T 



a o J |p(a)a 2 — k 



(3.175) 



Proof. Define t & = t(a = a©) then 



*© = / , _ _ + h, (3.176) 
\] \p( a ) a2 - k 



where to is the time now. The age of the universe is therefore given by 

r° da 

' a & ^\p{a)a? - k 



r° da , m . 

T = to - t Q = / . (3.177) 



□ 



Corollary 1. Assuming the universe starts with a big bang, the age of the universe r 
is the time elapsed since a = 0, and r is given by, 

r a ° da 

r= / - (3.178) 
^lp(a)a 2 - k 

Corollary 2. Assuming the universe starts with a bounce or a past sudden singularity, 
the age of the universe r is the time elapsed since a minimum value a = a*, and t is 
given by, 

f ao da 

r= / (3.179) 

Theorem 25. The age of the universe r is real and finite, that is, t < oo and r Gl, 
if only if the integral 

ao da . 

(3.180) 



is both real and convergent. 
That is, we demand both: 
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1. for reality: 

3k 

p(a) > — ; (3.181) 
cr 

2. for convergence: 

there exist constants a < 1, c > and a t > such that: Va G (a Q , a t ) 

3k 



and Va G (a e , a ) 



Proo/. 1. First, r G R if 



that is, if 



p(a) > ^ + cV Q - 2 ; (3.182) 
a z 



p(a) > p(a t ). (3.183) 



da 



a © J\p(a)a 2 - k 



(3.184) 



p(a)a 2 /3 - fc > 0, (3.185) 
which leads to the reality condition: 

p(a) > iL (3.186) 
a^ 

2. Second, r < oo if there exist constants a < 1, c > and a e > such that: 

Va G (a , a e ) 

y/p{a)a 2 /Z - k > ca a > 0, (3.187) 

that is 



p(a)>^ + cV a - 2 ; (3.188) 



a 



and, if furthermore Va G (a e , a ) 



Vp(a)a 2 /3 - k > ^p{a e )a 2 e /3 - k, (3.189) 

that is 

p(a) > p{a t ). (3.190) 

□ 
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3.4.2 Remaining lifetime of the universe T 

Definition 12. The remaining lifetime of the universe T is the remaining time from 
now, where a = oq, until a^ = a en d- 

Here a em j can represent a big rip or eternal expansion at a finite non zero rate 
(in this case a en d = oo), or a maximum value (for a turnaround or future sudden 
singularity or eternal expansion asymptotic to zero velocity in this Case ftend — ^max J • 



Theorem 26. The remaining lifetime of the universe T is given by, 

" a ° nd da 



I \p(a)a? — k 

Proof. Define t(aoo) the time when = a cm }, then 



(3.191) 



t{ aoo )= / _ + 1 . (3.192) 

|p(a)a 2 — k 

Therefore, the remaining lifetime of the universe T is given by 

*(aoo) - *o = / ^ i = (3.193) 
a o J±p(a)a 2 -k 

□ 

Corollary 3. Assuming the universe ends with a big rip, the remaining lifetime of the 
universe T is the remaining time from now a = a$ until a en< j = oo, and T is given by, 

da 



a <) J \p(a)a 2 - k 



(3.194) 



Corollary 4. Assuming the universe continues to expand at a finite nonzero rate for 
all eternity then, the remaining lifetime of the universe T diverges. Here T is the 
remaining time from now a = ao until a en d = oo, and T is given by, 

f°° da , 
T= / (3.195) 

\ \p(a)a 2 — k 
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Corollary 5. Assuming the universe ends by reaching a maximum value (possible 
turnaround or future sudden singularity), the remaining lifetime of the universe T is 
the remaining time from now a = oq until a en d = a max; ond T is given by, 

_ /" amax da 

T= - (3.196) 

J ao 



|p(a)a 2 — k 



Theorem 27. The remaining lifetime of the universe T is real and finite, that is, 
T < oo and T 6 1, i/ and only if the integral 

a cnd At. 

(3.197) 



J |p(a)a 2 — k 



is both real and convergent. 
That is, we demand both 

1. for reality: 



2. for convergence: 



p(a) > ^; (3.198) 



(a) case a en d = oo: If there exist constants (3 > 1, c > and a e > such that: 
Va G (a e , a cnd ) 

p(a) > + c 2 a 2/3 " 2 . (3.199) 
a 2 

(b) case a cnd = a max : If there exists constants a < 1, K > and a e > swc/i 
t/jat: Va G (a e , a max ) 

p(a) > ^ + 3(am g 2 ~ a)2 " > 4^ + 3(a -% a)2 ° . (3.200) 

a K a a max ^ a max 

And finally, Va G (ao, a e ) 

p(a) > p{a e ). (3.201) 
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Proof. 1. First, T G R if 



s J\p{a)a?-k 



(3.202) 



that is, if 

p(a)a 2 /3 -k> 0, (3.203) 
which leads to the reality condition: 

3A- 

p(a) > -j. (3.204) 



2. Second, 



[a) case a on d = oo: T < oo if there exist constants /3 > 1, c > and a e > 
such that: Va G (a e , a en d) 



y/p(a)a?/3-k > ca^ > 0, (3.205) 

that is 



3A; 



p(a) > ^ + c 2 a 2l3 ~ 2 . (3.206) 



a 



(b) case a en d = «max : T < oo if there exist constants a < 1, K > and a e > 
such that: Va G (a e , a max ) 

l/WI < 71 ^ _w ( 3 - 207 ) 



that is, 



> ^ + 3 ( a — fl ) 2Q > ^ + 3{a Z a 2 X 2 (3.208) 
a" 1 A^a 2 a* A^ai 



max max 



And, if furthermore Va G (ao, a £ 



y/p(a)a 2 /3 - k > V 'p(a e )a 2 t /3 - k, (3.209) 

that is 

p(a) > p{a e ). (3.210) 

□ 
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3.4.3 Total age of the universe T tota i 

Definition 13. The total age of the universe T total is the time elapsed since a = a Q 

Until doo = a e nd- 

Lemma 7. The total age of the universe T tota \ is the sum of the age of the universe r 
and the remaining lifetime of the universe T , that is, 



Lemma 8. The total age of the universe T tota i is real and finite if and only if: 

1. the present age of the universe is real and finite r£l and t < oo, 

2. and the remaining time of the universe is real and finite T £ R and T < oo. 

Theorem 28. The total age of the universe T to tai is given by, 



Theorem 29. The total age of the universe T tota i is real and finite, that is, T to tai < oo 
and Ttotai £ ^> if an d on ^V if the integral 



Ttotai = t + T. 





(3.212) 



Proof. Theorem l2"%l is deduced from Theorem |2U and Theorem EfD 



□ 




(3.213) 



is both real and convergent. 
That is, we demand both 



1. for reality: 




(3.214) 



2. for convergence: 
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(a) if there exist constants a < 1, c > and a e > such that: Va G (a , a e ) 

p( fl )>^ + cV tt - 2 ; (3.215) 
a z 

^ i. case a en d = oo: // there exist constants (3 > 1, c > and a e > sncn 
£/ia£: Va G (a e , a en d) 

p( a ) > ^ + C V^ 2 . (3.216) 

M. case a cn( j = a max : // inere exists constants a < 1, K > and a e > 
stzcn £/iat - Va G (a e , a max ) 

> ^ + 3 ( a "^ - fl ) 2a > ^ + 3(a "7 ; ar . (3.217) 

Proof. Theorem I2TH is deduced from Theorem 1231 and Theorem 1271 □ 

Now, applying the generic Theorem |2H1 to the five cases where the square root 
a/ p(a)a 2 /3 — k is real, we obtain the following corollaries: 

Corollary 6. Assume that the square root is real only on [0,a max ] and a max > clq. 
If there exist constants a < 1, c > and a e > such that: Va G (0, a e ) 



p{a) > — + c 2 a 2Q - 2 > — ; (3.218) 



3fc o 2rv — 2 

a z a 

and i/ taere exists constants a < 1, K > and a e > snca iaai: Va G (a e , a 



(a) > 3fc + 3(a a) 2 " > _3fc_ + 3(a a) 2 " . ^ 
a z A^a^ ai„„ A^ai 



"max "^max 



then 

Ttotal < OO ^ ^total G R. (3.220) 

Corollary IH1 deals with a universe that starts with a big bang (a Q = 0), and tends 
to a maximum size (a en d — Omax) m a finite time. This universe could possibly have a 
turnaround and contracting phase until a big crunch a = or endure a bounce a = a*. 
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Corollary 7. Assume that the square root is real only on [a*, oo]. 

Here a = is not possible, hence there is a bounce or a past sudden singularity. 
If there exist constants a < 1, c > and a e > such that: Va G (a*, a e ) 

p(a)>^ + cV«- 2 >^; (3.22 
ar a z 

and if there exists constants (3 > 1, c > and a e > swc/j t/iat: Va G (a e , oo) 



p(a)>^ + C V^ 2 >^; (3.222) 
cr cr 

Total < oo and T total G E. (3.223) 

Corollary [7| deals with a universe that starts with a bounce or a past sudden 
singularity (minimum value a & = a*), and ends with a big rip (infinite value a en d = oo) 
in a finite time. 

Corollary 8. Assume that the square root is real on the entire positive axis [0, oo]. 
If there exist constants a < 1, c > and a e > such that: Va G (0, a e ) 

p(a)>^ + cV a - 2 >^; (3.224) 
a z a z 

and if there exists constants [3 > 1, c > ana 1 a e > siica t/iat: Va G (a e , oo) 



p(a)>^ + c 2 a^- 2 >^; (3.225) 
a z a z 

then 

Ttotai < oo and T to tai G E. (3.226) 

Corollary |S] deals with a universe that starts with a big bang (a Q = 0) and ends 
with a big rip (a en d = oo) in a finite time. 

Corollary 9. Assume that the square root is real only on a finite interval bounded 
away from zero, [a*, a max ] and a* < a < a max . 

1/ t/iere exzsi constants a < 1, c > and a e > s«cn taat: Va G (a*, a e ) 

p(a)>^ + c 2 a 2 «- 2 >^; (3.227) 
a z a z 
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and if there exists constants a < 1, K > and a t > such that: Va G (a e , a maa: ) 

p(a) > ^ + 3(a — ~ a)2Q > 4^ + 3(a -% a)2 " ; (3.228) 

a A a a max a max 

then 

^totai < oo and T total G JR. (3.229) 

Corollary El deals with a universe that starts with a bounce or a past sudden 
singularity (minimum value a = a*) and ends with a turnaround or a future sud- 
den singularity (maximum value a en d = a max ). This universe could possibly have a 
turnaround and contracting phase until a big crunch a = or endure a bounce a = a*. 

Corollary 10. Disjoint unions of various of the above. 

This corollary corresponds to some of the cases described in the above corollaries. 

This case is not physically interesting since it corresponds to several unrelated 
universes that are totally independent of each other. 



3.4.4 Summary 

In this section, we have presented definitions and theorems regarding the conditions 
that should be placed on the the density as a function of the scale factor in order for 
the integral 




(3.230) 



to converge. We obtain model independent results that place constraints on p(a) under 
a minimum of technical assumptions. Whether or not the total age of the universe is 
a real and finite number depends on whether these conditions on the density p(a) hold 
or not. 
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3.5 Results and discussion 



In this chapter we have explored three issues concerning singularities in a FRW universe 
using only classical general relativity and an absolute minimum of technical input 
assumptions: 

• First we have developed an extensive catalogue of the various "cosmological mile- 
stones" found in the literature in terms of a generalized power series expansion 
of the FRW scale factor. 

If in the vicinity of any cosmological milestone, the input scale factor a(t) is a 
generalized power series, then all physical observables (H, q, the Riemann tensor, 
etc.) will likewise be generalized power series, with related indicial exponents 
that can be calculated from the indicial exponents of the scale factor. Whether 
or not the particular physical observable then diverges at the cosmological mile- 
stone is "simply" a matter of calculating its dominant indicial exponent in terms 
of those occuring in the scale factor. 

Indeed, for some unspecified generic cosmological milestone, that is defined in 
terms of the behaviour of the scale factor a(t), and which occurs at some finite 
time t Q , we can write in the vicinity of the milestone: 

a(t) = c \t - t |"° + ci]t - t @ \ m + c 2 \t - t |" 2 + c 3 \t - t^ + . . . (3.231) 

The rji are generically real and ordered such that t]q < rji < 772 < 773 . . . Also 
Co > though there are no a priori constraints on the signs of the other q, 
except by definition q ^ 0. 

This generalized power series expansion is sufficiently general to accommodate 
all commonly occurring models considered in the literature. Specifically, 

— ?7o > holds for the class of cosmological milestones such as big bangs or 
big crunches; 

— r]o < holds for the class of cosmological milestones such as big rips; 

— 7/0 = holds for the class of cosmological milestones such as sudden singu- 
larities or extremality events (bounces, turnarounds, inflexions); 
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Note that for extremality events the scale factor can be modeled using or- 
dinary Taylor series with rji G Z + . 



Second, with the notion of a generalized power series in hand, it is possible at 
a purely kinematic level to address the question of when a "cosmological mile- 
stone" corresponds to a curvature singularity, and what type of singularity is 
implied. 

In section 13. 2. 2\ we have shown that the Hubble parameter H(t) blows up in the 
vicinity of cosmological milestones such as big bangs/crunches and big rips (The- 
orem^, however it is finite for a specific class of sudden singularities/extremality 
events (Theorem HJ), specifically, 



lim H 

t-ft Q 



+00 


ci/c 

sign(ci)oo 
—00 



Vo > 
7?o = 

Vo = 
r/o = 

7/0 < 0. 



vi > i; 
m = !; 
771 e (0,1; 



(3.232) 



where £ Q is the the time of the event. 

In the same section, we have also demonstrated that the deceleration parameter 
q(t) is always finite for bangs, crunches and rips (770 7^ 0) but it is infinite only 
for a certain subset of the sudden singularities 



lim q = I 7/0 = 0; 771 = 1; 772 > 2; 

%e(i,2), 

(3.233) 

where £ Q is the the time of the event. 

We have considered the Riemann tensor and analysed its behaviour, in this 
way it is possible to classify all "cosmlogical milestones" as to whether they are 
polynomial curvature singularities or not: there are only a very few cases of cos- 
mological milestones that are not polynomial curvature singularities. In Theorem 



(1 - 77o)/r/o 






sign(ci[l - 771]) 00 


r/o = 







7/o = 


771 = 1 


— 2C2Q)/ c\ 


7/o = 


771 = 1 


-sign(c 2 )oo 


r/o = 


771 = 1 
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El we mention that those cases correspond to sudden singularities/extremality 
events (t]q = 0) and asymptotically Milne big bangs/crunches. In detail, the only 
cases where cosmological milestones are not polynomial curvature singularities 
are: 

— Vo = 771 > 2 or; 

— Vo = ??i = l V2 > 2 or; 

— r)o — 1 771 > 3 and fc = — 1, Co = 1. 

However, when looking at time derivatives of the curvature tensor, most of these 
limited cases disappear: the only two situations in which a cosmological mile- 
stone is not a derivative curvature singularity are if: 

— 7/0 = 0, r\i G Z + ; corresponding to an extremality event; 

— 7/0 = 1, = —1, c = 1, rji G Z + , and 7?! > 3; corresponding to a FRW 
geometry that smoothly asymptotes near the cosmological milestone to the 
Riemann-flat Milne universe. 

• Third, this definition of cosmological milestones in terms of generalized power 
series enables us to perform a complete model-independent check on the validity 
or otherwise of the classical energy conditions. 

In particular we provide a complete catalogue of those bangs/crunches, sudden 
singularities and extremality events for which the NEC, the WEC, the SEC and 
the DEC are satisfied. 

Depending on one's attitude towards the energy conditions [3], one could use 
this catalogue as a guide towards deciding on potentially interesting scenarios to 
investigate. In particular, 

— The NEC, the WEC and the SEC hold independent of the sign of space 
curvature k all the way down to the singularity, for big bangs/crunches for 
which < 7/0 < 1. That is, in the vicinity of the singularity the dominant 
term in the scale factor a(t) is bounded from both above and below by 




(3.234) 
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— "Big rips" (770 < 0) will violate the NEC in the vicinity of the milestone 
and therefore, they will also violate every other energy conditions. 

— The SEC will be violated for all "violent" bangs and crunches with 770 > 1 , 
though the NEC and the WEC will still be satisfied in this case. 

— The DEC holds independent of the sign of space curvature k all the way 
down to the singularity, for big bangs/crunches for which 1/3 < 770 < 1. 
That is, in the vicinity of the singularity the dominant term in the scale 
factor a{t) is bounded from both above and below by 

co(t - Q < a dominant (t) < c (t - i*) 1/3 . (3.235) 

— The DEC is also satisfied for a limited range of sudden singularities (7/0 = 0). 
Sudden singularities of this type are sufficiently "gentle" that at least for 
some of them even the DEC can be satisfied all the way to the singularity: 
a conclusion completely in agreement with Barrow and Tsagas jS]. 

Writing the scale factor as a generalized power series expansion has allowed us to 
derive important and very general results on the various cosmological milestones both 
on a kinematical and dynamical levels. These theorems have been developed with an 
absolute minimum of technical imput. Simply by calculating the dominant indicial 
exponent of the scale factor we have been able to determine whether or not some 
particular physical observable (expressed as power series) diverges at the cosmological 
milestone. 
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Chapter 4 
Gravastars 



This chapter introduces the reader to the idea of gravastars, how they are configured, 
as well as some of their properties. We present results deduced from these properties 
but we do so from an agnostic point of view as to the existence or non-existence of 
gravastars. 

The study and understanding of the features of gravastars are important to com- 
prehend what observational data regarding astrophysical black holes are telling us. 

4.1 Definitions 

4.1.1 GRAvitational VAcuum STARS 

The concept of black hole 1 is generally accepted in the general relativity community, 
but one sometimes encounter some scepticism concerning the reality of the mathemat- 
ical solution and wariness regarding the interpretation of observational data jT]. The 
Schwarzschild solution presents a central singularity at r = and an event horizon at 
the Schwarzschild radius 2 RschwarzschUd = 2M. 

The GRAvitational VAcuum STAR (gravastar) model is a proposal by Emil Mot- 
tola and Pawel Mazur EH EH| to replace black holes. Instead of having a star 

1 See section 12.5.11 for more details. 
2 See section 12.5.21 for more details. 
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contracting and collapsing until matter arrives at a singularity at the centre, the 
gravastar model suggests that a gravitationally collapsing star would force spacetime 
itself to undergo a phase transition that would prevent further collapse. Thus the star 
would be transformed into a spherical "quantum vacuum" surrounded by a form of 
super-dense matter. 

It has been speculated that violent creation of a gravastar might be an alternate 
explanation for gamma ray bursts. Gravastars could also be a solution for the black 
hole information paradox 3 . The Gravastar is theorized to have very low amounts of 
entropy, on the contrary, a black hole apparently has a billion times more entropy 
than the star it formed from. A problem with the theory regarding the creation of a 
gravastar is whether or not a star would be capable of shedding enough entropy upon 
implosion. Externally, a gravastar appears similar to a black hole: it is visible only 
by the high-energy emissions it creates while consuming matter. Astronomers observe 
the sky for X-rays emitted by infalling matter to detect black holes, and a gravastar 
would produce a similar signature. 

4.1.2 Common models 

In the Mazur-Mottola model, the concept of Bose-Einstein condensation 4 in gravita- 
tional systems is extended to compact object with an interior de Sitter space (with 
an equation of state p = —p > 0) , and an outer region of the gravastar consisting 
of a (relatively thin) finite-thickness shell of stiff matter (p = p). The exterior of the 
gravastar is surrounded by a Schwarzschild vacuum geometry (p = p = 0). 

In addition to these three layers, the model requires two infinitesimally-thin shells 

3 The black hole information paradox results from the combination of quantum mechanics and 
general relativity. In 1975, Stephen Hawking showed that black holes should slowly radiate away 
energy, which poses a problem. From the no hair theorem one would expect the Hawking radiation 
to be completely independent of the material entering the black hole. However, if the material entering 
the black hole were a pure quantum state, the transformation of that state into the mixed state of 
Hawking radiation would destroy information about the original quantum state. This violates the 
rules of standard quantum mechanics and presents a seeming physical paradox. 

4 In Bose-Einstein condensate, all matter (protons, neutrons, electrons, etc..) goes into what is 
called a quantum state creating a "super-atom" . 
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with surface densities cr±, and surface tensions $±. These compensate the discon- 
tinuities in the pressure profile and stabilize this 5-layer construction, introducing 
delta-function anisotropic pressures [S3 EHl EH] • 

4.1.3 Problem 

Is it possible to replace the thin shell (from the Mazur-Mottola model) completely 
with a continuous layer of finite thickness? For physical reasons, we find it useful to 
minimize the use of thin shells as they really are a mathematical abstraction. 

Is the pressure anisotropy (implicit in the Mazur-Mottola infinitesimally thin shell) 
a necessity for any gravastar-like objects? Is it really possible to build a gravastar-like 
objects using only perfect fluid with a continuous layer of finite thickness? Could a 
horizon or naked singularity form and why? 

Assuming pressure is continuous and differentiable, we will first analyze the result- 
ing static geometry with the resulting isotropy TOV equation (with isotropic pressure) 
and then the resulting anisotropy TOV equation (anisotropic pressure). 

4.2 Properties 

4.2.1 A static spherically symmetric geometry 

To describe the geometry of the gravastar, we use spherical coordinates (t, r, 6, 0) and 
we assume that the geometry is static, and spherically symmetric. 

• A static solution means that there exists a hypersurface orthogonal Killing vector 
that is timelike near spatial infinity. This implies that the spacetime metric ds 2 
can be chosen to be invariant under a time reversal about any origin of time, 
e.g. all cross terms dx t dx l V i vanish in the chosen coordinate system. 

• A spherically symmetric solution means that there exists a privileged point, e.g. 
the origin, such that the system is invariant under spatial rotations about the 
origin. 
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The most general static spherically symmetric line element in four dimensions can 
be written in the canonical form |79j : 

ds 2 = -f(r)dt 2 + h(r)dr 2 + r 2 dQ, (4.1) 

where dVt 2 = d6 2 + sin 2 6d(j) 2 . 

For any spherically symmetric stress-energy tensor, we have: 

Tab = pu a u b + p r n a n b + p t (g ab + u a u b - n a n b ) , (4.2) 

where p r is the radial pressure, pt the tangential pressure, u a the fluid 4-velocity 
pointing in the same direction as the static Killing vector field £ a (in order to be 
compatible with the static symmetry spacetime) and n a is a unit radius vector. This 
property of u a leads to a relation between the 4-velocity and the function / appearing 
in the general metric: 

u a = -f 1 2(dt) a (4.3) 

Now, introducing the Einstein's equations developed in section I2~3*l we obtain relations 
between elements of the stress energy tensor and the functions / and h from the metric. 
Specifically, the equation in the tt components, using "geometrized units" 5 and with 
orthonormal components, involve: 

G U = 8vrT £t - = Sirp (4.4) 
Equation f)4.4j) can be rewritten involving the function h: 

~(r(l-h- 1 ))=8np. (4.5) 

Hence, the solution for h can be calculated: 

. . / 2m(r) \ 1 . 
h(r) =1 K -L , (4.6) 



where 6 

pr 

m(r) = Air p{f)f 2 df. (4-7) 



o 



5 Geometrized units mean that G and c are set equal to one, see section ETUI for more details. 
6 Note that in order to avoid a "conical singularity" in the metric at r = 0, (in other words if 
we require the geometry to be regular at r = 0), then the integration constant mo in the equation 



m(r) = 47r J£ p(f)r 2 dr + mo must be set equal to 0. 
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Therefore, the tt-field equation yields equation (J4.7|) . The term m(r) represents the 
total mass-energy inside the radius r integrated from the centre r = until arbitrary 
r. Note that if r = cannot be defined (e.g., in the case of wormhole throats), the 
integration would be from an arbitrary constant tq to r. Also note that this equation 
of the mass is precisely the equation you would expect for a spherically symmetric 
mass in an Euclidean space. 

The rr-field equation from the relation ()2.15|) of Einstein field equations, in ge- 
ometrized units and with orthonormal components, gives: 

Grf = BnTff = 87ip r . (4.8) 

Equation (|4.8j) involves the gu = —f{f) term of the metric, but it is more physically 
convenient to express f(r) as a function of the quantity g(r) defined below. 

Definition 14. Let g(r) be related to the gu component of the metric as follows: 

9(r) = \j-H-9u) (4-9) 
= a r g„ (4.10) 

where a r is the radial acceleration of an observer "at rest" , which means that the 4 
velocity vector is parallel to the Killing field 7 . 

Now, after rearranging this definition, we can write the g tt component of the metric 
in terms of g(r) : 

gu = -exp\2^ g(r)drj . (4.11) 
Now, after further calculations, equation (|4.8|) . that is the G rr equation implies that 

r 2 i _ Mr) x 



The physical meaning of g(r) is that this function represents the locally measured grav- 
itational acceleration, which is pointing inwards for positive g(r). Note that equation 
(J4.7|) also means that 

^ = 4vrpr 2 . (4.13) 
dr 



7 In a Killing vector field, (r, 9, <j>) is constant. 
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Definition 15. We define the average density p = -j — -, so that the gravitational 



acceleration g(r) can be written: 



4?r 3 
3 ' 



g(r) = f^^. (4,4) 

r 

Finally, adopting (t, r, 8, <fi) coordinates we write any static spherically symmetric 
geometry, and in particular for gravastars, in the form 8 

/ r 00 \ dr 2 
ds 2 = -expl2 g{f)dr\dt 2 + - ^— + rW*. (4.15) 

v 1 ~) 

Of course, the metric of the exterior of the gravastar is obviously the Schwarzschild 
exterior solution (|2.31|) 9 . 



4.2.2 A TOV equation 

The remaining field equation Ggg from the Einstein field equations (|2.3|) . gives an 
expression for the radial pressure p r , however, the rather messy algebra can be cir- 
cumvented by the use of the Bianchi identity ()2.9|) to replace it with the covariant 
conservation equation of the stress-energy tensor: 

V b T ab = 0. (4.16) 

The stress-energy tensor is given in equation ()4.2j) . after differentiating, this yields to 

- r = -{p + Pr)g + • (4.17) 

dr r 

Equation f)4. 17|) is known as the anisotropic Tolman-Oppenheimer-Volkoff (TOV) 
equation. In the case of isotropic pressure (p = p r = pt) this leads to the more familiar 
isotropic TOV equation: 

% = -(P+P)9, (4-18) 
8 This spacetime geometry is also referred to as an "interior solution" . 

9 Rcnicmber that the unique vacuum solution to the Einstein's equations is the Schwarzschild 
exterior solution by Birkhoff stated in Theorem ^ 
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which can be rewritten, using Definition an d the expression of the average density 
p, as 

dp _ (p + p)(m + Airpr 3 ) _ An r (p + p)(p + 3p) (4 19) 

dr r 2 [1 — 2m(r)/r] 3 1 - 2m(r)/r 

Definition 16. Let's define the dimensionless anisotropy parameter A for anisotropic 
pressures by 

A = (4.20) 
P 

Therefore, in terms of this dimensionless parameter, the anisotropic TOV equation 
becomes! 

dp r _ Anr (p + p r ) (p + 3p r ) | 2pA 

<ir 3 1 — 2m{r)/r r 

The TOV equation shows that pressure contributes as a source of gravity, and one can 
notice that it follows that the pressure is higher than in Newtonian stars of the same 
density profile. It also means that for a given p(r), the central pressure p c required for 
equilibrium is always higher in general relativity than in Newtonian theory. Therefore, 
it is harder to maintain equilibrium in general relativity. 

How to solve the TOV equation? 

For full determination of the variables p r , pt, p, m, one needs the TOV equation ()4.21|) 
but also 

• equation (j4.13|) which gives the derivative of the mass dm/dr as a function of 
the density p, 

• an equation of state for the radial pressure of the stellar material 

p = p(p r ) (4.22) 

• an equation of state for the tangential pressure of the stellar material 

P = P(pt) (4-23) 

• boundary conditions: 
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— spacetime must be everywhere locally Lorentzian and in particular at the 
centre of the star, it is therefore necessary that the mass vanishes at the 
centre 10 : 

m = at r = (4.24) 

— at the star's surface the interior spacetime geometry (|4.2.1|) must join 
smoothly to the exterior Schwarzschild geometry (|2.31)l . therefore 

m = M at r = R sur face, ( 4 -25) 
and by definition, the radial pressure p r is zero at the surface of the star. 
The method to produce a model is straightforward providing that: 

• the equations of state mentioned before are specified, 

• central density p c or central pressure p c 11 for the star are given. 

One can integrate the coupled hydrostatic equilibrium equation (|4.21|) and the mass 
equation (j4.13|) outward from the centre, beginning with the initial conditions m = 
and p r = p c at the center. The integration terminates when the pressure falls to zero 12 
at the surface of the star. The value of the mass at this surface radius represents the 
total mass-energy M that appear in the Schwarzschild exterior solution (|2.31j) . 

Qualitative differences introduced by A. 

Let's have a qualitative look at the difference that anisotropy introduces in the TOV 
equation with dimensionless anisotropy parameter A. Remember that if we assume 
the WEC 13 holds, then p ^ everywhere in the star. 

10 This condition (locally Lorcntz) is actually responsible for leading exactly to equation l|4.7fl . that 
is m(r) = Jg Airr 2 pdf. 

11 Spherical symmetry requires that p r = p t = p c at the centre. 

12 Note that for gravastars, theoretically, the pressure starts off negative, passes through a first zero, 
increases to a maximum pressure before falling down to zero at the surface. 
13 See section EOTD 
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1. case A > 0: (e.g. p r < p t ) 

Typically, this condition holds for rigid solids, in telluric planets like the earth 
for example. If the radial pressure is negative, and equivalently pushing inwards, 
the tangential pressure, created by the rigidity of the surface, can be responsible 
for preventing the object from collapsing. In other words, the tangential pressure 
is doing more to stop gravitational collapse than the radial pressure. 

2. case A < 0, and < p t < p r : 

This particular case occurs in fluids with weak surface tension and where gravity 
plays an important role to hold the fluid together. More importantly, the radial 
pressure must be greater at the centre than it would be in a completely isotropic 
object as the pressure profile is steeper. 

3. special case: p t < 0: 

Objects satisfying this condition are fluids with strong surface tension, such as 
raindrops for example. The radial pressure must also be greater at the centre 
than it would be in a completely isotropic object as the pressure profile is steeper. 
However, in this case, the radial pressure does not necessarily go to zero smoothly 
at the surface. 

4.2.3 Specific key features 

In this section, we present the class of spacetime geometries that we are particularly 
interested in for the gravastar model. 

In the same spirit of Mazur and Mottola (23 EH EH]> and Laughlin et al |2lH l21j. 

we assume that the weak energy condition holds through out the configuration (i.e. the 
density is positive), but we permit the pressure to become negative in the gravastar 
interior. The density does not need to be continuous 14 but to avoid infinitesimally 
thin shells (purpose of this model), one must demand that the radial pressure p r is 

14 Typically, the density is not continuous at the surface of the gravastar. 
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continuous. The radial pressure is drawn qualitatively in Figure 14.11 




Figure 4.1: Qualitative sketch of radial pressure as a function of r for a gravastar . 

To "smooth out" the infinitesimally thin shells of the Mazur and Mottola gravastar 
model, we consider static spherically symmetric geometries such that: 

• Inside the gravastar, r < R, the density is everywhere positive and finite. 

• The central pressure is negative, p c < 0, and in fact p c = —p c . 
(We do not demand p = —p r = —p t except at the centre 15 .) 

Positive energy density but negative pressure, and origin free of any mass sin- 
gularity, are the characteristic features of de Sitter space. The idea of replacing 
the Schwarzschild singularity with de Sitter vacuum goes back at least to 1965. 
Gliner |37] interpreted p = —p as corresponding to a "vacuum" and suggested 
that it could be a final state in a gravitational collapse. 

15 Rcmcmber that because of the spherical symmetry p r — pt at the centre. 
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Indeed for the the specific metric 



ds 2 = -c 2 7 (r)rft 2 + 7 - 1 (r)dr 2 + rW, 



(4.26) 



and the choice of cosmological-constant-like matter leads to 




(4.27) 



where r s 



2GM 



with M the total mass of the object, and with 




) 



2 



r < r 



S ! 



7(r) 



(4.28) 



r > r 



S ) 



Note that the interior de Sitter region may be also interpreted as a cosmological 
spacetime, with the horizon of the expanding universe replaced by a quantum 
phase interface. Also note the unpleasantness at r = r s : Laughlin et al call this 
a "phase transition" ; Mazur-Mottola seek to smooth out the unpleasant region. 

Hence, the characteristics of the gravastar at the centre, are those of a de-Sitter 
spacetime geometry. 

• The spacetime is assumed to not possess an event horizon. 
This implies that Vr we have 2m(r) < r. 

These three features, positive density, negative central pressure, and the absence of 
horizons, are the three most important features characterizing a gravastar. 
Other important features are: 

• To keep the centre of the spacetime regular, we enforce both p',(0) = and 
Pc = p r (0) = Pt(0). 

• There should be a pressure maximum in the general vicinity of the Schwarzschild 
radius 16 , r max ps -Rschwarzschiid, satisfying p r (r m3X ) > 0, and p^(r max ) = 0. 

(This permits the physics in the region r ^> r max to be more or less standard.) 
16 See section 12X21 
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• There should be exactly two radii where the radial pressure vanishes: 

— The first pressure zero p r (r ) = 0, where p' r (r ) > 0, and 

— the second pressure zero p r (R) = 0, where p' r (R) < 0. 

The point R, (which by construction must satisfy R > r max > r ), is called 
the surface of the gravastar. 

• The pressure profile p r (r) should be continuous. 

(In contrast, it is sometimes useful to allow pt{r) to be discontinuous.) 

• The strong energy condition [SEC; p + p r + 2p t > 0] is definitely violated, at 
least near the centre of the gravastar. 

At the centre of the star, p c = p r = p t = —p with p > 0, therefore 

p + Pr + 2p t = p- p-2p= -2p< 0. (4.29) 

• We choose to enforce the null energy condition [NEC; p + Pi > 0] throughout the 
gravastar. In view of our first comment that density is everywhere positive, this 
implies that we are enforcing the weak energy condition [WEC; p + Pi > and 
p>0]. 

• We impose no restriction regarding the dominant energy condition [DEC; p > 
and \pi\ < p], and in fact we shall see that the DEC must fail in parts of the 
gravastar that are sufficiently "close" to forming a horizon. 

Why does the radial pressure goes through two zeros? 

It is clear on Figure l4~T1 that our model describes a radial pressure that starts negative 
at the centre, then increases till p r = 0, becomes positive, and starts decreasing after 
a maximum p r = p max at r = r max eventually reaching a second zero at the surface. 

By definition, when the pressure is zero it corresponds to the surface of the star 
so we could have smoothly joined the vacuum interior to an exterior Schwarzschild 
solution (|2.31J) . The reason for not doing so is a purely pragmatic one based on the 
fate of infalling positive-pressure matter. If we start without a positive-pressure region 
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of type (vq,R), then any infalling positive-pressure matter that accumulates above r 
will automatically generate a positive-pressure region of type (r ,R). The only way 
to avoid a positive-pressure region of type (r ,R) is if the negative-pressure matter 
in the (0, ro) region immediately catalyzes any infalling positive-pressure matter into 
negative-pressure matter. In this entire chapter, we will keep the negative-pressure 
matter deep in the core, discretely hidden behind a layer of positive-pressure matter. 

How does our model relate to others? 

We have highlighted specific features that a perfect fluid sphere gravastar, with contin- 
uous pressure, should have. But how does our gravastar model relate to other features 
in previous models? 

• In the Mazur-Mottola model , the infinitesimally thin shell interface is the lim- 
iting case where r max — > tq, while p = —p r = —pt is strictly enforced for r < r^. 
Also, an additional thin shell is placed at the surface so that the radial pressure 
is positive on the left limit of the surface radius p r (R~) > 0. 

• A simplified variant of the Mazur-Mottola 3-region model is also considered in 
|77j but the surface radius tends to r max , e.g. R — > r max . So in this model 
?"o = r max = R arid there is a single thin shell at r max with de Sitter geometry 
inside and Schwarzschild geometry outside. 

• In the Laughlin et al model, r -> i?sch W ar ZSC hiid horn below, while R -> -Rschwarzschiid 
from above. And in the region r < -Rschwarzschiid > the equality p = —p r = —p t 
remains all the way through. There is implicitly a singular infinitesimally thin 
shell located exactly at -Rschwarzschiid with infinite surface tension. 

• The Gliner jSH] and Dymnikova PIES EE1 EH EB] proposals all satisfy the 
constraint p = —p r everywhere throughout the configuration. We assume p = 
—p r only at the centre, this equation does not need to hold everywhere. 
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4.3 The problem of the negative-pressure perfect 
fluid sphere gravastar. 



In this section, we will only consider perfect fluids, and therefore the pressure is now 
isotropic p = p r — p t , which means that the dimensionless anisotropy parameter is 
null throughout this whole section A = 0. The TOV equation is described by equation 
(|4.19jl in this section. 



4.3.1 Four inconsistencies 

In this section, we will have a look at four cases where some inconsistency within the 
isotropic TOV equation arise for perfect fluid sphere gravastars. 

First inconsistency 

Theorem 30. The isotropic TOV equation cannot hold at the point tq where the 
radial pressure is zero p{tq) = 0. 

Proof. By assumption, the LHS of the TOV isotropic equation, dp/dr, is positive. As 
shown on Figure (|4.1j) . the pressure is increasing in the region r < r . However, we 
can calculate the exact value of the pressure gradient at r = ro where p = 0: 



dp 
dr 



47rrn op , 

(4.30) 



= ro 3 l-2m/r ' 

The density is positive (as we assume the WEC) and therefore the average density p 
is positive as well. Furthermore r < -RschwarzscMd, which implies that 2m/r < 1. As 
a consequence, it follows that the RHS of equation ()4.30|) is negative. 

Therefore, there is an inconsistency, the LHS cannot be positive while the RHS is 
negative. Thus, the isotropic TOV cannot hold at the point r . □ 

Second inconsistency 

Theorem 31. The isotropic TOV equation cannot hold at the point r max where the 
radial pressure is positive and reaches a maximum p(r max ) > 0. 
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Proof. Now let's consider the point of maximum radial pressure r max . We can calculate 
the value of the TOV equation at this point: 



dp 
dr 



47rr max (p + p)(p + 3p) 



2m/r E 



At r max the pressure is a maximum and positive for the class of models we consider, 
therefore p(r max ) > and its derivative is zero dp/dr\ r=r = 0. But the density is 
also positive everywhere p > by assumption, therefore the LHS is zero while the 
RHS is negative. 

This is another inconsistency, the isotropic TOV cannot hold at the point r max . □ 
Third inconsistency 

Theorem 32. The isotropic TOV equation cannot hold in the region r < r < r max 
where the radial pressure is positive and increasing. 

Proof. Now if we look at the region between the first and second inconsistency, r < 
r < r max; we can notice that the object under investigation has an increasing positive 
radial pressure, therefore the LHS of the TOV equation is positive in this entire region. 
However, because we assume p > which also implies that p > 0, and that there is no 
horizon for at least r < r max , the RHS of the TOV equation is negative for this entire 
region. 

_4 ! rr jp + pKp + Sp) < Q 
3 1 — 2m jr 

This is the third inconsistency and it follows that isotropic pressure in that interval 
r o ^ r ^ r max is not able to satisfy the TOV equation, and therefore, we conclude that 
a static spacetime geometry can only be obtained with the introduction of tangential 
pressures, that is, anisotropic pressures. □ 

Fourth inconsistency 

Definition 17. Let's define r g the location where the quantity p + 3p in the isotropic 
TOV equation \4.19\) changes sign. 
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Theorem 33. The isotropic TOV equation cannot hold in the region r g < r < r 
where the radial pressure is negative and increasing and with (p + 3p) rg - 

Proof. The same arguments hold for a larger interval below ro and into the negative 
pressure region. 

Assuming the WEC through the entire configuration implies that NEC 17 holds, 
from which follows that 

p + p^O Vr:r<r . (4.33) 
But remember that the LHS of the TOV equation is positive, therefore 

^k±lM±M <0 =, ? + 3 p<0 Vr<r„ (4.34) 
3 1 — 2m jr 

Indeed, at the centre of the gravastar, we have 

(p + 3p) c = p c + 3 Pc = -2 Pc < 0. (4.35) 

However, exactly at the point ro, where the pressure is zero, we have 

(p + 3p) ro =p ro >0. (4.36) 

Therefore, it means that the term p + 3p changes sign somewhere in the interval 
< r < r . 

In the region r g < r < ro, we have both: 

p + 3p > and ^ (4.37) 
dr 

This is another inconsistency. Therefore, the pressure isotropy fails for the entire 

region r g < r < ro and quite possibly fails for an even larger region. □ 

Qualitative definitions 

The following definitions relate the qualitative physics to the different regions of the 
gravastar described in the previous sections. 



17 The WEC implies the NEC, see sections EOTD and ETO 
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Definition 18. The region < r < r g , where the physics is qualitatively similar to 
that of de Sitter space, will be referred as the "core" . In the core, the local acceleration 
of equation due to gravity is outward. 

Definition 19. The region r g < r < r max , where the physics is still definitely "un- 
usual", will be referred as the "crust". In the crust, the local acceleration due to 
gravity is inward, but the pressure still rises as one moves outward. 

Definition 20. The region r max < r < R, where the physics is still definitely "nor- 
mal", will be referred as the "atmosphere" . In the atmosphere, the local acceleration 
due to gravity is inward, and the pressure decreases as one moves outward. 

Figure (|4.2j) presents a qualitative sketch of the previous defintions. 

With these definitions, we see that pressure is guaranteed to be anisotropic through- 
out the "crust" . Note that even if we were to dispense with the entire positive-pressure 
region by chopping the gravastar off at r , there is still an anisotropic crust in the re- 
gion (r g ,r ]. 

In conclusion, we have shown that there are inconsistencies within the isotropic 
TOV equation from which follows that a static spherically symmetric object with 
positive density, negative central pressure and vanishing pressure at the surface cannot 
be supported by isotropic pressures alone. 

There are no perfect fluid gravastars. 

4.3.2 A particular gravastar case study: the Schwarzschild 
interior solution and the Buchdahl-Bondi bound 

The Schwarzschild interior solution illustrates what can go wrong when trying to build 
a perfect fluid gravastar. 

The Schwarzschild solution in standard coordinates (t, r, 9, <fi) is given by equation 
(|2.31|) . But the general Schwarzschild interior solution for values of r smaller than the 
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r = r g r r max R 

Figure 4.2: Qualitative sketch of gravastar labelling the "core", "crust", and "atmo- 
sphere" . 

Schwarzschild radius -Rschwarzschiid is given by: 



2 



if = -{ A + B^l- 2 -^f> j *» + ^ T -i 3Er I *" + rW (4.38) 

where dfl 2 = d9 2 + sin 2 9d(p 2 and m*(r) = M(r/R) 3 . This general solution describes 
static spherically symmetric stars with constant density independent of r from the 
centre to the surface of the star. 

Using boundary conditions, we can determine A and B as functions of the surface 
radius R and the total mass M of the star. The density can be calculated with the G a 
term of the Einstein tensor and gu part of the metric (|A.7|) due to the field equations 
flZED: 

8ir P *g tt = G tt (4.39) 
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The value of the density is constant and depends only on M and R that are fixed for 
a specific star. 

The radial pressure can also be calculated exactly with the field equations in the rr 
components as a function of A and B from the general metric. However, the expression 
of p* is quite complicated and messy, but the important feature is, that the surface 
radius R can be deduced as a function of A and B from the boundary condition: 

at Surface p*(R)=0 => f x {A,B,R, M) = 0. (4.41) 

Furthermore, the total mass of the star is: 

rr=R(A,B,R,M) 

M = 4vr / p*r 2 dr f 2 {A,B,R,M) = 0. (4.42) 

Jr=0 

Therefore, solving the following system 1)4.43)1 of two equations for A and B 

f 1 (A,B,R,M) = 
f 2 {A,B,R,M) = 

leads to A(R, M) and B(R, M) which determines the metric of the Schwarzschild 
interior solution entirely as a function of the total mass M and the surface radius R. 
Indeed, it can be found that: 



2M 

A = \ 1 



B = 

In this case the metric f)A.7|) becomes: 



R (4.44) 




^-h/l-Y-^ 1 -^! ^+(-^t)^ 2 + ^ 2 M5) 

where m*(r) = M(r/R) 3 . 

In 1916, Schwarzschild integrated exactly by hand the isotropic equation ()4.19)) 
yeilding 

3y/l - 2M/R -y/1- 2m,(r)/r 
where M and R are the total mass and surface radius. Note that the denominator of 
the pressure is exactly 3gu from the tt-component of the metric (j4.45j) . 
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Central pressure p* and the Buchdahl-Bondi bound 

The central pressure required for equilibrium of a uniform density star is 

1 - y/j - 2M/R 

V c — P* — ; (4.47) 

3v/l-2Af/i?- 1 

For i? > M, equation ()4.47jl reduces to the Newtonian value ()4.49|1 18 19 : 

P := (4.49) 

The sign of the central pressure p* is determined by specific conditions only on the 
denominator as the numerator is always positive: 

2M 



1 - y/1 - 2M/R > ^> > — (always true). (4.50) 

R 

The central pressure p* c in the gravastar becomes infinite when 

2M 



Zy/\-2M/R = \ =► -H-=g. ( 4 - 51 ) 



and therefore, p* > when 



2M 



3 V / 1-2M/ J R>1 ^ ~^ < g- (452) 

This result on the compactness, \ = 2M/ R, is known as the Buchdahl-Bondi bound: 
the compactness of a static spherically symmetric fluid x is bounded above by 8/9. 
This is quite significant as it shows that a star cannot get arbitrarily close to forming 
a black hole (% = 1). 

However, in our particular gravastar model, we deal with negative central pressures 
p* < 0, and this implies that we have the relation: 

8 2M 

p* c <0 <=► -<_<!. (4.53) 



IS 



Note that the isotropic TOV equation in a Newtonian star is 



dp / \ / \ m(r) , . 4 . , 

■£ = -g(r)p(r) = ~^p(r) = -Trrp. (4.48) 



19 Value given in [Sj 
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Of course, while this is a perfectly sensible solution in the mathematical sense one 
would generally rule it out physically because of the negative central pressure. If we 
leave our prejudice against negative pressure aside, adopting the gravastar philosophy, 
we will find some interesting results coming from the choice of (I4.53J) . 



First order pole 

Theorem 34. For a Schwarzschild interior solution, with constant positive density 
p* > 0, and a negative central pressure p* < 0, the pressure profile given by equation 
\4 ■ 4 6| ) will have a first order pole at 



8/9 

= 3R\I 1 - ^M/R and f urthermore r p°ic < R - ( 4 - 54 ) 

Note that Theorem EI] assumes constant positive density p* > and a central 
negative pressure p* < 0, but it does not assume that at the centre p* = —p*. As a 
matter of fact, the only way to obtain exactly p* = —p* at the centre is if R — 2M 
that is if there is a Schwarzschild radius. Indeed, p* = —p* if and only if: 



p* 1 - Jl - 2M/R 

-1, (4.55) 



that is, if, 



finally, if 



P* 3y/l-2M/R-l 

2y/l - 2M/R = 0, (4.56) 
R = 2M. (4.57) 



Proof. From the pressure profile given in equation (|4.46|l . there is a pole when its 
denominator cancels out or, in this case, when the tt-component of the metric (|4.45jl 
is zero. 



3^/1-2^-^1-27^,(^ = 0, (4.58) 
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where m*(r) = M(r/R) 3 . We can extract the radius r from this relation and obtain a 
value for r at the pole. 



Mr 2 M , . 

Mr 2 / M 



1_2 1P~ = 9 V 1_2 liJ < 460 > 

Mr 2 / M\ , 

-2—— = 8-9x2— (4.61) 



R 3 \ R) 



= 3R ^-mii (463) 

Now we can see what the conditions to have the upper bound R on r po i are: 



3R ^~wk < R < 464 > 



>/»-5^g) < o (*■«) 

2M 

==>- — — < 1 true for no event horizon (4.67) 
R 

Therefore, for a Schwarzschild interior solution, with negative central pressure and 
constant positive density, that has no event horizon, there exists a pole such as 

r P oic < R, (4.68) 

where R is the surface. □ 

Note that as 2M/R goes from 8/9 to 1, the position of this pole moves from the 
centre of the star to the surface of the star. The situation is qualitatively sketched in 
figure 14. M 
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Figure 4.3: Qualitative sketch of pressure pole in the interior Schwarzschild solution 
for 2M/R > 8/9. 

Theorem 35. The NEC must be violated sufficiently close to the pole at r po \ e = 

QD A 8/9 
6n \l 1 2M/R ■ 

Proof. The compactness is bounded 8/9 < 2M/R < 1 which implies that the denomi- 
nator term in equation (j4.46J) is negative, and therefore, that the pressure p* tends to 
— oo as r approaches r polo . 

The NEC holds, in the vicinity of the pole, if and only if: 

p*+p*>0. (4.69) 

However, near the pole, 

lim (p* +p*) = -oo. (4.70) 

' ^pole 

Therefore, the NEC must be violated sufficiently close to the pole. □ 
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Curvature singularity 

We are going to see now that Schwarzschild interior solution gravastars are clearly 
unphysical, even if one is willing to accept negative pressures and even violations of 
the NEC. It is unphysical because the pressure pole implies a curvature singularity. 
With the metric given by equation ([4.45)1 . we can calculate the Riemann tensor in an 
orthonormal basis R^^. In particular, for isotropic pressures, we know that, 

Rt* = -^ = ^(P + 3 P ) (4.71) 
Qttgee o 

Now, we know that those specific terms of the Riemann tensor will have a problem 
at r = r po i c . As shown previously for the pressure pole, the term g tt is zero at r = r po i c 
and 

g tt = 0{[r-r polc ] 2 ), (4.74) 
whereas g rr is finite at r = r pole and 

g rr (r polc ) = 9 (1 _ \ MjR) > 0. (4-75) 
Concerning the trir-term of the Riemann tensor, we have 

lim R trtr = 0, (4.76) 

'* ^' pole 

dBsjr = 3 2_M\ { ~ 2Wr) 
r-r polc dr 81R 3 V R J (1 - 2M/R) K 1 

Hence, we have for the trtr-term of the Riemann tensor at the pole: 

lim Rtrtr = O {[r - r polc \) . (4.78) 

^* ^polc 

Therefore, the tftf-term of the Riemann tensor becomes infinite at r = r po i e , as 
gu becomes the dominant term: 

lim R m = oo. (4.79) 

^* ^polc 
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On a similar way, we have gee and g&s finite at the pole, and 



9ee(r po ie) = 9R 2 [l- -f^- ) > 0, (4.80) 



2M/R 

^(r po ie) = 9 sin 2 6R 2 (l - J > 0. (4.81 

Concerning the t9t9 and t(j)t(j) terms of the Riemann tensor, we have, 



lim Rtete 


= o, 


^* ^pole 




, . dRtete 
lim — - — 


3 


^»-poi c dr 


R 


lim Rt<f>t<t> 


= o, 


T ^polc 




lim 


3 


->r pole dr 


R 



(4.82) 



(4.84) 

2 / a /o \ V2 



2M/RJ 



sin 2 (5) . (4.85) 



Hence, we have for the and terms of the Riemann tensor at the pole: 

lim Rtete = O ( [r - r polc ] ) , (4.86) 
lim R Ht<t> = O ([r - r pole ]) . (4.87) 

^* ^polc 

Therefore, the i§i9-term and the t0t0-term of the Riemann tensor become infinite at 
r = r P oie, as g a becomes the dominant term: 

lim Rieie = oo» ( 4 - 88 ) 

^polc 

lim = a)- (4-89) 

T ^polc 

However, the rest of the non-zero orthonormal components of the Riemann tensor are 
finite, 

RfM = 2 ~§ ( 4 - 9 °) 
R*M = ^ ( 4 - 91 ) 

2M 

Rem = #r ( 4 - 92 ) 
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To conclude, we have shown that 3 orthonormal components of the Riemann tensor 
(Rifir Q4.79|) . Rf§t§ (|4.88|L R§i§i (|4.89|) ) are infinite at r = r po i e , and therefore, the 
pressure pole implies a naked singularity. 

The reason why we mention this specific example is because we shall soon see that 
this "pressure pole" behaviour is generic. Continuous solutions with isotropic pressure 
are not possible 20 . 



4.3.3 Arbitrary gravastar case study: the fate of a negative 
pressure perfect fluid sphere 

In this section, we consider an arbitrary gravastar as defined in section (|4.2|) . and we 
will show that if the gravastar configuration is perfect fluid, the same problems as 
those encountered in the previous section ()4.3.2j) will arise in this more general case. 

For any arbitrary gravastar that is static spherically symmetric, remember that 
the geometry is given by 



/ r°° \ dr 2 
exp ( 2 / g{r)dr\ dt 2 + + r 2 dQ 2 , (4.93) 

v 1 ~) 

= m(r)+4ir Pr (r)r3 



ds* 

where 



Theorem 36. If an arbitrary gravastar is perfect fluid and finite, and if we permit 
the violation of the NEC, then there exists a pressure pole of order 1 and it is the only 
way the TOV equation can be satisfied. 

Proof. Using equation ()4.13|) that relates the density p to the mass m(r), we can obtain 
an expression for the compactness \ = 2m{r)/r which satisfies: 

" - P] = ^ [3(p +p)-(p+ 3p)]. (4.95) 





2m(r) 


' 2m 8- 




= 87r p r - — = — 




r 





20 Consider for example ^0], where a perfect fluid Chaplygin gas p oc 1/p is considered. The surface 
of their configuration occurs at p — where p = — oo, at least as one approaches the surface from 
below. So the surface of their configuration is a naked singularity, in agreement with the observations 
above. 



116 



But the first term on the RHS is non-negative by the NEC, while the second term is 
by definition negative on [0, r g ), so the compactness 2m(r)/r is monotone increasing 
on the range [0, r g ). This is somewhat unusual, while the compactness of a normal 
perfect fluid star tends to increase as one moves outwards, in a normal star it also can 
be subject to oscillations that make the overall picture quite subtle |4*2| 178] . 

Let's consider now a perfect fluid sphere with negative central pressure that satisfies 
the NEC. Since we have already seen that isotropy is violated on (r g , r max ), the only way 
we can maintain the perfect fluid nature of the sphere is if r g — > oo, which also implies 
that ro — > oo. But since we do not want a horizon to form, the compactness must be 
bounded above by unity. But we have just shown the compactness is monotonic and 
therefore, 

Um 2m(r) = ^ x ,£(0,l]. (4.96) 

r— too T 

So not only does a NEC-satisfying perfect-fluid gravastar expand to infinite volume, 
it also has infinite mass. 

To avoid the physical size of the gravastar blowing up to r — > oo, our options now 
are rather limited: we could permit the development of a horizon at finite r, which 
defeats the whole point of the exercise, or we could permit something even worse. 

First order pressure pole 

If we permit NEC violations then it is possible to arrange for the development of a 
pressure pole at finite r < r . Let's see how this is possible by first considering the 
following definition. 

Definition 21. Set 

p(r) « (4.97) 
r — r p 

in the vicinity of the pressure pole of order 1 at r p where T is a positive constant. 

This definition is compatible with the isotropic TOV equation. First for r < r p the 
pressure is negative, and for r > r p the pressure is positive, which is appropriate for a 



117 



modified gravastar model. Secondly 



V « 7 ^ < 0. (4.98) 



Third, assuming p remains finite, close to the pole (r m r p ) we have 

(P + p)( m + 47rpr 3 ) [p + r/(r — r p )] [m + 4nr 3 T/ (r — r p )] 

r 2 (l — 2m/r) r 2 (l — 2m(r)/r) 

[T/(r-r p )] [47rr 3 r/(r - r p )] 
r 2 (l — 2m(r)/r) 



(4.99) 
(4.100) 



47rr T 2 1 

<0. (4.101) 



(1 - 2m(r p )/r p ) (r - r p ) 2 
So the TOV equation can be satisfied in the vicinity of the pole provided 

(p + P) ( m + 47rpr 3 ' 1 



r 2 (l — 2m/r) 

that is provided we set 



pj. at r = r g , (4.102) 



1 - 2m(r p )/r p > Q 



4irr p 



That is, if the gravastar configuration is perfect fluid and finite in extent, then this 
pole is the only way the TOV equation can be satisfied. □ 



n 



th order pressure pole 



Now we can look at higher order pole and see if they also have the nice property of 
being compatible with the isotropic TOV equation. 

Definition 22. Set 



P(r) « 7 yi, (4-104) 

(r - r p ) 

in the vicinity of the pressure pole of order n > 1 at r p where T is a positive constant. 



This definition is compatible with the isotropic TOV equation for odd value of n 
only. First for r < r p the pressure is negative if n is odd, and for r > r p the pressure 
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is positive if n is odd. However, this definition is not compatible with the isotropic 
TOV equation for even values of n. Secondly 

— TlY 

P' ~ 7 u^TT < °- ( 4 - 105 ) 

Third, assuming p remains finite, close to the pole (r « r p ) we have 

(P +p)(m + 47rpr 3 ) [p + Y / (r — r p ) n ] [m + Anr^Y / (r — r p ) n ] 



r 2 (l — 2m/r) r 2 (l — 2m(r)/r) 

[17 (r - r p ) n ] [4vrr 3 r/ (r - r p ) n ] 



r 2 (l — 2m(r)/r) 

A-KTpY 2 1 



(4.106) 
(4.107) 



(1 - 2m(r p )/r p ) (r - r p ) 2n 
So the TOV equation can be satisfied in the vicinity of the pole provided 

(p + P) { m + 47rpr 3 ' 1 



< 0. (4.108) 



r 2 (l — 2m/r) 

that is provided we set 



V r at r = r g , (4.109) 



1 - 2m(r p )/r p 



4irr p 



n (r -r p ) n_i ^ 0. (4.110) 



First, equation ()4.110j) is in total contradiction with the Definition E21 of Y that is 
supposed to be a constant. Secondly, it shows that the isotropic TOV equation cannot 
hold at the pressure pole of order n. 

Therefore, higher-order poles do not even have this nice property of being compat- 
ible with the isotropic TOV equation. 



Theorem 37. An arbitrary finite gravastar with a perfect fluid configuration that has 
a pressure pole leads to a curvature singularity. 

Proof. For any isotropic static spherically symmetric spacetime the orthonormal com- 
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ponents of the Riemann tensor are (see, for instance, [741 p 110]): 





An r . 


fp)-2p]; 


(4.111) 


Rieie = 


An , 

T |3p + 


pi; 


(4.112) 


Rr6r6 = 




•Pi; 


(4.113) 


R§$94> = 


f M- 




(4.114) 



It is clear that if the pressure has a pole of order n ^ 1, then, providing p remains 
finite as per assumption, as Rfcfc and depend on p, they will also have a pole, 
and consequently, a naked singularity will arise. □ 

Even if the isotropic TOV equation holds with pressure pole in the gravastar con- 
figuration, this situation is physically inappropriate as the pressure pole introduces a 
naked singularity. The physically correct deduction form this analysis is that gravastar- 
like objects must violate pressure isotropy. 

Summary 

In this section, we have shown that a static spherically symmetric object with positive 
density, negative central pressure, and vanishing pressure at the surface, (which are 
the defining features of a gravastar), cannot be supported by isotropic pressure alone 
without a pole in the pressure. Despite the fact that the presence of a simple pole in the 
pressure is compatible with the isotropic TOV equation, we must reiterate that such a 
pressure pole is unphysical because it is a naked singularity. Therefore, gravastar-like 
objects must violate pressure isotropy. 

4.4 Gravastars with anisotropy 

In section (14. 3 j) . we have demonstrated that gravastars cannot have isotropic pressure 
alone, and that, as a consequence, anisotropic pressure must be part of the config- 
uration or at least in some parts of it. In this section, we are going to see when 
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the anisotropy for the pressure is necessary for the TOV equation to hold, and what 
specific features the equation of state must have. 



4.4.1 Bounds on the anisotropy pressure 

Once we accept that perfect fluid spheres are not what we are looking for to model 
gravastars, one might wonder what happens to the Buchdahl-Bondi bound for isotropic 
fluid spheres. It has been shown that for p' < and pt < p r the 8/9 bound still holds 21 . 
However if the transverse stress is allowed to exceed the radial stress, (e.g., p t > p r ), 
then the upper limit shifts to 2M/R < k < 1, where k depends on the magnitude of 
the maximal stress anisotropy r 



Remember that in section EQ] we defined the crust of the gravastar as the region 
r g < r < r max , where the physics is still definitely "unusual", and where the local 
acceleration due to gravity is inward, but the pressure still rises as one moves out- 
ward. The dimensionless anisotropy parameter A is a function of the radial pressure, 
tangential pressure and positive density: 



Pt ~ Pr r 



p \ p 



(4.115) 



p 2 

Theorem 38. In the crust of a gravastar, the compactness is not limited to the 
Buchdahl-Bondi bound, but only by the magnitude of the maximal pressure anisotropy 
and the regularity of the metric: 

2m / V 

A > '—— > r„ < r < r (4.116) 

~ 4 1 - 2m/r 9 v ' 

and 

rp' . rpL 1 2m /r 

0<-^<A<-^ + '—— r < r < r max 4.117 

~ 2p 2p 4 \-2mjr K ' 

Furthermore, the tangential pressure is greater than the radial pressure p t > p r for 

r g <r < r max . 



21 See section (|Qlfl) 

22 For more related work on anisotropic stars see |45l 1441 
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Proof. In section R.3. 11 we have determined that anisotropy is necessary in the crust 23 , 
that is, when r g < r < r max . We can now find bounds on A in that interval by inserting 
the definition of the gravitational acceleration g(r) into equation (|4.115|) . We get, 



Pr 
P 



p r \ m + 47rp r r 
P 



^ _ 2m 



(4.118) 



Now we know that for the interval r £ [?"o, r max ], the radial pressure is positive 
p r > and increasing so p' > 0. Therefore, we can write the inequality, 



A * 2 ' 

A > 



P, 



m 



p J r 2 [l — 2m/r] 



> 



m 



2 r 2 [l - 2m/r] 
that leads to the simple lower bound: 



> 0, 



(4.119) 
(4.120) 



2m /r 
4 1 — 2m jr 



> 0. 



(4.121) 



The inequality on A (J4.121|) . for the interval r G [r , r max ], shows that it is positive, 
consequently, in this region, we have pt > p r > 0. 

Now let's see what happens for the interval r £ (r g , ro). We know that in this area, 
the radial pressure is negative p r < 0, but it is increasing and dp r /dr > 0. From these 
conditions on the radial pressure and its derivative, we can write, 



A ^ 
A ^ 
A ^ 



rip; 



P 



m + 4npr r 3 
r 2 (1 — 2m/r) 
m 



rp r r 

~2p~ + 2 r 2 [l - 2m/r] 
1 2m jr 



rp' r 



+ 



(4.122) 
(4.123) 
(4.124) 



2p 4 1 - 2m/r ' 

Remember that we assume that the NEC holds in the configuration of the gravastar. 



2:s 



Isotropy possibly fails for an even larger region. 
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This leads to 



p + p r > 



Pr > 
Pr 

P 



> 



-P 
-1 



1 + 



Pi 



> 0. 



(4.125) 
(4.126) 
(4.127) 

(4.128) 



This last inequality (|4.128J) . allows us to write a lower bound on A: 

m + Anpr r 3 



^ + 1 



A>^>0. 



p J r 2 (l — 2m/r) " 2p 

Consequently, for the interval r £ (r g , ro), we find the weaker bounds 

o < LEz < a < rp ' r 1 1 2m//r 



(4.129) 



2p 



2p 4 1 - 2m/r 



Theorem 39. In the interval [ro,r max ] ; if the DEC holds, then 

DEC A < 1. 



(4.130) 
□ 

(4.131) 



Proof. For the first inequality on A (|4.121jl . in the region r £ [ro,r max ], the DEC is 
satisfied if 



p^O, 

"P < Pr ^ P, 
-p^Pt^ P- 



This leads to the more useful inequalities, 

Pr 
P 

Pt 
P 



^ 1, 



C 1. 



(4.132) 
(4.133) 
(4.134) 



(4.135) 
(4.136) 
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But because the radial pressure is positive for the interval [ro,r max ], we can write, 

P P 

And therefore, if the DEC is to be satisfied, it implies the following upper bound on 
A: 

DEC A < 1. (4.138) 

□ 

Theorem 40. Any gravastar that is sufficiently close to forming a horizon, that is 
2m/r > 4/5, will violate the DEC in its "crust". 

Proof. If the DEC is to be satisfied we must at the very least have A < 1. Therefore, 
when the DEC holds, we have 

1 > 

1 > 

4 ^ 
4 
5 

Hence, the DEC is guaranteed to be violated whenever 2m/r > 4/5. If the gravastar is 
sufficiently close to forming a horizon, in the sense that 2m jr > 4/5 somewhere in the 
range [ro, r max ], then the DEC must also be violated at this point. Even if we were to 
discard the entire positive-pressure region (tq,R), we can nevertheless still apply this 
bound at r itself: if 2m(r )/r > 4/5 then the DEC is violated at r . Consequently, 
any gravastar that is sufficiently close to forming a horizon will violate the DEC in its 
"crust" . □ 

4.4.2 Minimizing the anisotropy 

We have shown that the anisotropy is necessary for gr avast ar-like objects in their 
"crust", where the pressure is increasing as one moves outwards and the local force 



1 2m / r 

A > '—— 

4 1 — 2m/r 

1 2m/r 



> 



4 1 — 2m jr 

2m 

5 — 

r 

2m 



> 



(4.139) 
(4.140) 

(4.141) 
(4.142) 
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of gravity is inwards. In this section, we attempt to minimize the region over which 
anisotropy is present. 

Remember that he dimensionless anisotropy parameter A is deduced from the 
anisotropy TOV equation (j4.21j) . and gives, 



A = Pt-Pr = r 



(4.143) 



What happens at r 9 ? 

At the point r g , we can write (since g(r g ) = by definition II 7|) 

A(r ff ) = r ^ rg) > 0. (4.144) 
To confine the anisotropy to the smallest interval possible, we have to set 

p' r (r 9 ) = 0, (4.145) 
which corresponds to an inflexion point for the radial pressure. 

What happens at r ? 

At the point r , the radial pressure is zero, and we can write, 

Hence, at the point of zero radial pressure, the anisotropy cannot vanish. 
What happens at r max ? 

At the point r max , the radial pressure reaches its maximum value. If we take the limit 
from below, we can write, 

A(w) = i , 2m J r - (l + P AU+ > • (4-147) 



4 1 - 2m/r max \ p J \ m 
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Remember that the anisotropy is not necessary beyond the peak V r > r max (in the 
atmosphere), and therefore, it is possible to arrange A = 0. To confine the anisotropy 
to the smallest interval possible, we have to set 

A(r-r+J = 0. (4.148) 

However, this leads to a discontinuity in the derivative of the radial pressure p' r and 
the anisotropy parameter A, as well as a "kink" in the pressure profile p r at r max . 

From the TOV equation (j4.21|) . we can deduce the limit from above for the pressure 
derivative, 



A at r~ ax can be rewritten as 

Therefore, we can still keep the radial pressure and the density continuous which will 
imply that, 

K(r max ) = A(r max ). (4.151) 

'max 

The implications of confining the pressure anisotropy to the smallest interval possi- 
ble are shown in figure PO] where the anisotropy is confined to the region r G (r g , r max ]. 



How does this relate to other gravastar models? 

The Mazur-Mottola model is recovered as the limiting case where r g — > tq <— r max 24 . 
This results in the fact that all the important anisotropy is confined in their inner thin 
shell. The anisotropy A — > oo, because p' — > oo. Effectively A is replaced by choosing 
an appropriate finite (in this case negative) surface tension •& and surface energy density 
a which is given by the Israel-Lanczos-Sen junction conditions |4*| IHT | IHT^ | 177]. 

The second, outer thin shell which is present in the Mazur-Mottola model is not 
a physical necessity, but is merely a convenient way to avoid an infinitely diffuse 
atmosphere that would arise otherwise from the equation of state p = p. A finite 
24 See section 114. 1.211 for more details 
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Figure 4.4: Qualitative sketch of radial pressure as a function of r for a gravastar that 
is as near as possible to a perfect fluid. Note the inflexion point at r g and the kink at 

''max' 

surface radius R can be modelled by altering the equation of state slightly to include 
a finite surface density ps which is reached for vanishing pressure: p(p) = p$ + p. 
Then, the outer thin shell can be omitted when joining the gravastar metric onto the 
Schwarzschild exterior metric. 



4.4.3 Features of the anisotropic equation of state 

In the case of a perfect fluid, the geometry is completely defined by the set of differ- 
ential equations (|4.1Hjl . (|4.19jl plus initial conditions and an equation of state, which 
is commonly written as p = p(p) or equivalently p = p(p). 

However, in this chapter we have pointed out that gravastars with perfect fluid 
sphere are a lost cause, and one needs to replace the isotropic TOV equation (j4.19J) 
with its anisotropic counterpart (|4.21j) . By doing so, we also need to introduce an 
additional free function A(r). In order to close the set of equations, it is necessary to 
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define two equations of state. 



Two equations of state 

Two equations of state must be chosen, but which criteria does one take into account? 

• p r (p) and p t (p) 

This is a rather strong assumption that forces p r and p t to change in lock-step. 

• density profile p(r) and a single equation of state / pressure profile p r (r) and a 
single equation of state 

This is the strategy adopted, for example, in [23J ESI EH EH EH] • 

• Specify any two profiles p(r), p r (r), and pt{r) by hand, and use the anisotropic 
TOV to calculate the remaining profile. 

This is rather unphysical. 

• finding an additional (differential) equation 

This is another possibility to obtain a well defined solution, that might be moti- 
vated by some appropriate variational principle. An example would consists of 
minimizing the "total anisotropy" 



No matter whether it is an equation of state for the transverse pressure or a 
variational principle for the total anisotropy, the extra equation should (if the 
gravastar model is to be even qualitatively correct) be responsible for maintaining 
stability of the gravastar over a wide range of total masses and central pressures. 
In other words, the closed set of equations must be self-regulating if it is to be 
physically interesting. The gravastar should shift to a new stable configuration 
when the total mass changes. 




(4.152) 
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The fact that, there will be regions where A > in the gravastar configuration 
suggests that it might be most efficient to choose the two distinct "equations of state" 
as being equations for p and for A. 



What variables should these equations of state depend on? It seems obvious that 
they should depend on the radial pressure p r , but in view of the lower bounds on A in 
equation (j4.116j) . it is clear that to avoid the formation of an horizon in the gravastar 
configuration, the equations of state should be sensitive to the compactness 2m/ r. 

Therefore, we should posit equations of state of the form 

p — p(r,p r ,2m/r); A = A(r,p r , 2m/r); (4.153) 

and solve the paired differential equations 

dp r (r) [p(r,p r (r),2m(r)/r) -\- p r (r)] [2m(r)/r + 8ir p r (r) r 2 } 

dr 2r [1 — 2m(r)/r] 

| 2 p(r,p r (r),2m(r)/r) A(r,p r (r), 2m(r)/r) _ 

r 

dm(r) 2 

— = 4-7T p(r,p r (r),2m(r)/r) r . (4.155) 

dr 

The success of this gravastar model depends on whether or not one can find phys- 
ically realistic equations of state that take into account the effect of the compactness 
2m /r for large ranges of central negative pressure p c and total mass M. 

Potential issue: 

Some traditional relativists might argue that there is an important issue of principle, 
by stating that the variable 2m /r is not detectable by local physics. There might also 
argue that an equation of state that depends on 2m /r somehow violates the Einstein 
equivalence principle 25 . 

But this is incorrect because 2m /r is certainly measurable by quasi-local effects 
in small but finite size regions. For any static spherically symmetric spacetime the 



25 See section El for the equivalence principle 
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orthonormal components of the Riemann tensor are (see, for instance, [73J p 110]): 

[3(p-p r + 2p t )-2p]; (4.156) 

[Spr + p]; (4.157) 

[3p-p]; (4.158) 

[2p\. (4.159) 

The Riemann tensor is certainly measurable in finite-sized regions, so in particular 
p is measurable. Likewise r is measurable in finite-sized regions, and therefore 2m/ r = 
(87r/3) p r 2 is measurable. 

Consequently the compactness 2m/ r is a quasi-local measurable quantity and it 
can meaningfully be put into the equation of state without violating the equivalence 
principle. 

Note that in either the Mazur-Mottola scenario EU EH] or the Laughlin et al. 
scenario [201 |2T] the gravastar material is assumed to be a quantum condensate, and 
therefore sensitive to non-local physics. The point is that we do not need to appeal to 
quantum non-locality to get 2m/ r into the equation of state. 

However, note that while this argument demonstrates that possible equations of 
state are at least conceivable, that is not the same as explicitly demonstrating that 
such equations of state actually exist. 

4.5 Results and discussion 

This chapter has hopefully served as a brief introduction to the idea of gravastars and 
how they are configured, as well as some of their properties. 

We have discussed the qualitative features one would expect from a gravastar con- 
figuration that has finite continuous radial pressure profile everywhere and that does 
not have delta function transition layers (|4.2|) . More importantly, we posit that 
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• Inside the gravastar, r < R, the density is everywhere positive and finite. 

• The central pressure is negative, p c < 0, and in fact p c = — p c , and the spacetime 
in the centre is regular. 

• The spacetime is assumed to not possess an event horizon. 

• We choose to enforce the WEC throughout the gravastar but impose no restric- 
tions on other energy conditions. 

• Our model describes a radial pressure that starts negative at the centre, then 
increases till p r = 0, becomes positive, and starts decreasing after a maximum 
Pr = Pmax at r = r max eventually reaching a second zero at the surface. 

Adopting (t, r, 9, 0) coordinates, the interior geometry for gravastars can be written 
in the form, 

/ f°° \ dr 2 

ds 2 = - exp 2 / g(r)df dt 2 + + r 2 dtt 2 (4.160) 

\ Jr J (1 - 2m(r)/r) 

where dfl 2 = d6 2 + sin 2 9d(p 2 , and m(r) represents the total mass-energy inside the ra- 
dius r integrated from the centre r = until arbitrary r. The gravitational acceleration 
g(r) is written as: 

9{r) = — l-2m(r)/r - (4161) 
We have used the isotropy TOV equation where the tangential pressure is equal to 
the radial pressure p r = p t , 

f r = -{p + P)9, (4-162) 
and the more general anisotropy TOV equation, 

dp r _ Attt (p + p r ) (p + 3p r ) | 2pA (4 163) 

dr 3 1 — 2m{r)/r r ' 

where we have defined the anisotropy dimensionless parameter 



Pt ~Pr 
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(4.164) 



From theses specific features and the isotropy TOV equation, we have analyzed and 
deduced that gravastars with perfect fluid sphere do not exist, that is, gravastar-like 
objects must have anisotropic pressure. 



We also have extracted some generic information concerning the features that the 
equations of state of a gravastar should have. More specifically, any gravastar-like 
object, that has negative central pressure p c , must have anisotropic pressure in their 
"crust"; the "crust" being the region where the pressure is increasing as one moves 
outward and the local force of gravity is inward 1)4.3.1)1 . Since the compactness is 
monotonic, we have, assuming the NEC, 

lim— = X.i X*e(0,l]. (4.165) 

r— >oo T 

Not only do perfect fluid gravastars lead to inconsistencies when the isotropy TOV 
equation holds for finite continuous pressure, but, when the pressure is allowed to be 
discontinuous for the TOV equation to hold, they lead to an infinite-size infinite-mass 
object, or a naked singularity as the pressure exhibits a pole 1)4.3.21 and 14.3.31) . 

If we assume the DEC, we can explicitly bound the magnitude of the anisotropy, 
required in the crust, in terms of the local compactness 2m(r)/r. We have shown that, 

A> r Q <r<r (4.166) 

~ 4 1 - 2m/r 9 u v ; 

and 

0< — <A<— + '—— r <r<r max . (4.167) 

~ 2p 2p 4 1 - 2m/ r V ' 

The second inequality (j4.167j) is weaker than the first one (|4.166|) . and it might still 
be possible to find better bounds. 

The magnitude of the anisotropy becomes arbitrary large for gravastars that are 
sufficiently close to forming a horizon, that is, 

2m 4 . 
— > -, 4.168 

r 5 

in which case the DEC must be violated f)4.4.1j) . 
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The compactness bounds are harder to state as well: if one is prepared to ignore 
the DEC then it seems that the compactness can be arbitrarily close to one. 

When trying to minimize the anisotropy region in the gravastar model, one has to 
pay the price of a discontinuity in p' r and A, and a "kink" in the radial pressure profile 

Pr- 



We have also deduced, that for configurations that are horizon-avoiding for any 
total mass and negative central pressure, the equations of state of gravastar matter 
must depend on the local compactness 2m(r)/r. We posit that the equations of state 
are of the form 

p = p(r,p r ,2m/r); A = A(r, p r , 2m/r); (4.169) 

and that one should solve the paired differential equations 

dp r (r) \p(r,p r (r),2m(r)/r) +p r (r)] [2m(r)/r + 8tt p r (r) r 2 ] 

dr 2r [1 — 2m(r)/r] 

| 2 p(r,p r (r),2m(r)/r) A(r,p r (r), 2m(r)/r) _ 

r 

dm(r) A _ f __ ^ o ^ t ! ^ 2 



dr 



47T p(r,p r (r),2m(r)/r) r . (4.171) 



The success of this gravastar model depends on whether or not one can find phys- 
ically realistic equations of state that take into account the effect of the compactness 
2m j V. However, we have not explicitly demonstrated that such equations of state 
actually exists. Note that, as discussed in ([4.4.3)) . this unusual equation of state does 
not violate the equivalence principle. 

Even though we are personally agnostic to the existence or non-existence of gravas- 
tars, the study and understanding of their properties are important to comprehend 
what observational data regarding astrophysical black holes are telling us. 
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Chapter 5 
Conclusion 



Two main problems in general relativity have been addressed in this thesis: 

• The first problem is relevant to cosmology and singularities. 

• The second looks at the gravastar model developed by Mazur and Mottola, a 
possible alternative to black holes. 

Main results have already been summarized to a large extent at the end of each 
respective chapter. However, the main important results are highlighted below: 

• Cosmology: 

Three issues have been explored concerning this topic: 

— We have developed a complete catalogue of the various cosmological mile- 
stones in terms of a generalized power series expansion of the FRW scale 
factor. This power series expansion is sufficiently general to accommodate 
all commonly occurring models considered in the literature. 

— With the notion of a generalized power series in hand, it is possible at a 
purely kinematic level to address the question of when a cosmological mile- 
stone corresponds to a curvature singularity, and what type of singularity 
is implied. An important result is that there are only a few cases of cosmo- 
logical milestones that are not polynomial curvature singularities. And the 
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only two situations in which a cosmological milestone is not a derivative 
curvature singularity correspond to : 

* an extremality event; 

* a FRW geometry that smoothly asymptotes near the cosmological mile- 
stone to the Riemann-flat Milne universe. 

— Finally, this definition of cosmological milestones in terms of generalized 
power series enables us to perform a complete model-independent check on 
the validity or otherwise of the classical energy conditions. In particular we 
provide a complete catalogue of those bangs/crunches, sudden singularities 
and extremality events for which the NEC, the WEC, the SEC and the DEC 
are satisfied. Depending on one's attitude towards the energy conditions, 
one could use this catalog guide towards deciding on potentially 
interesting scenarios to investigate. In particular, the DEC is satisfied all 
the way to the singularity for a few range of sudden singularities (r]o = 0) 
that are sufficiently gentle. 

• The gr avast ar model: 

— We have discussed the qualitative features one would expect from a gravas- 
tar configuration that has finite continuous radial pressure profile every- 
where and that does not have delta function transition layers. 

— From these specific features and the isotropic TOV equation, we have an- 
alyzed and deduced that gravastars with perfect fluid sphere do not exist, 
that is, gravastar-like objects must have anisotropic pressure. 

— We also have extracted some generic information concerning the features 
that the equations of state of a gravastar should have: for configurations 
that are horizon-avoiding for any total mass and negative central pressure, 
the equations of state of gravastar matter must depend on the local com- 
pactness 2m(r)/r. The success of this gravastar model depends on whether 
or not one can find physically realistic equations of state that take into 
account the effect of the compactness 2m(r)/r. 
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To conclude, I would like to point out several open problems that have not yet 
been analyzed. 

The results posited in cosmology assume an isotropic and homogeneous FRW uni- 
verse. However, it would be interesting to look at what happens when assuming an 
inhomogeneous universe. Calculations would become much harder and not as simple 
as in the homogeneous case, but will we still have some energy conditions satisfied all 
the way through the cosmological milestones? 

Concerning the gravastar model, its success depends on whether or not one can 
find physically realistic equations of state that take into account the effect of the 
compactness. However, we have not explicitly demonstrated that such equations of 
state actually exist. This is a challenge that we leave for the future. 
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Appendix A 

Main spacetime metrics considered 



This appendix assembles all the main metrics that have been mentioned in this thesis. 

A.l Static spacetime 

For a static spacetime, in the chosen coordinate system with arbitrary {x a }, the metric 
components are of the form: 

3 

ds 2 = -V 2 (x\ x 2 , x 3 )dt 2 + 9ab(x\x 2 , x 3 )dx a dx\ (A.l) 

a,b=l 

where V 2 = — £ a £ a and £ a is a Killing vector field. 
This metric is first mentioned in section 12.5.21 

A.2 Schwarzschild exterior solution 

The Schwarzschild solution, describing static spherical symmetric vacuum spacetimes, 
is in standard coordinates (t,r,9,4>): 




(A.2) 
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where dVL 2 = <i# 2 + sin 2 9d<p 2 . This solution is also known as the "Schwarzschild exterior 
solution". This is true for any spherically symmetric vacuum solution to Einstein's 
equations. M is a parameter that can be interpreted as the conventional Newtonian 
mass that would be measured by studying orbits at large distances from the gravitating 
source. 

This metric is first mentioned in section 12. 5. 21 



A. 3 Friedmann- Robertson- Walker geometry 

In cosmology, the homogeneous and isotropic cosmological model called the Friedmann- 
Robertson- Walker geometry, (FLRW Friedmann-Lemaitre-Robertson- Walker geome- 
try) is a good description of our universe and is derived by symmetry considerations 
without using the Einstein equations of general relativity. The metric is given by: 

ds 2 = -dt 2 + a(t) 2 | + r 2 [d6 2 + sin 2 9 dtf] j ( A.3) 

where, a(t) is the scale factor of the universe. There are only three values of interest 
for the parameter k: 

• k = — 1, this corresponds to a negative curvature (for the hyperboloid) 

• k — 0, this corresponds to no curvature (flat space) 

• k — +1, this corresponds to a positive curvature (for the 3-sphere) 

The assumptions of homogeneity and isotropy alone determine the spacetime metric up 
to three discrete possibilities of spatial geometry k and the arbitrary positive function 
of the scale factor a(t). 

This metric is first mentioned in section ETT1 
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A. 4 Static spherically symmetric spacetime 



The most general static spherically symmetric line element in four dimensions can be 
written in the canonical form |79| : 

ds 2 = -f{r)dt 2 + h{r)dr 2 + r 2 dQ, (A.4) 

where dQ 2 = dQ 2 + sin 2 9d(f> 2 . 

This metric is first mentioned in section 14.21 



A.5 Interior of a gravastar 

Adopting (t, r, 9, 0) coordinates we write any static spherically symmetric geometry, 
and in particular for the interior of a gravastar, in the form: 

/ r°° \ dr 2 

ds 2 = - exp 2 / g(r)df dt 2 + -. r + r 2 d£l 2 , (A.5) 

V I / M - 2 n^l\ 

where the gravitational acceleration g(r) is: 

4vrrp + 3p r (r) 
9( r ) = — - 2 m(r) • ( A - 6 ) 

r 

This general solution is also know as an "interior solution". Of course, the metric of 
the exterior of the gravastar is obviously the Schwarzschild exterior solution (jA.2|h 
This metric is first mentioned in section 14.21 



A. 6 General Schwarzschild interior solution 

The Schwarzschild exterior solution in standard coordinates (t, r, 8, <fi) is given by equa- 
tion (jA.2|) . A different solution is the general Schwarzschild interior solution for values 
of r smaller than the Schwarzschild radius -Rschwarzschiid is given by: 

ds 2 = - [ A + B Jl- ^M^j dt 2 + { \-^\ dr 2 + r 2 dQ 2 (A.7) 
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where dVL 2 = d9 2 + sin 2 6d<f) 2 and m*(r) = M(r/R) 3 . This general solution describes 
static spherically symmetric stars with constant density independent of r from the 
centre to the surface of the star. 
This metric is first mentioned in section 14.3.21 

A. 7 Schwarzschild interior solution as a function of 
the total mass M and the surface radius R 

If one wishes to determine the metric of the Schwarzschild interior solution entirely as 
a function of the total mass M and the surface radius R instead of A and B, it can 
be found that: 




(A.8) 



In this case the metric becomes: 





) 



2 



) 



ds 2 



1 



2M 



2m* (r) 



dt 2 + 



1 



dr 2 + r 2 dn 2 (A.9) 



R 



r 



1 



2m, (r) 



r 



where m*(r) = M(r/R) 3 . 



This metric is first mentioned in section 14.3.21 
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Appendix B 



Necessary and sufficient conditions 
for big bangs, bounces, crunches, 
rips, sudden singularities, 
extremality events, and more... 

Celine Cattoen and Matt Visser 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

Until recently, the physically relevant singularities occurring in FRW cosmologies had tradi- 
tionally been thought to be limited to the "big bang" , and possibly a "big crunch" . However, 
over the last few years, the zoo of cosmological singularities considered in the literature has 
become considerably more extensive, with "big rips" and "sudden singularities" added to 
the mix, as well as renewed interest in non-singular cosmological events such as "bounces" 
and "turnarounds". In this article we present an extensive catalogue of such cosmological 
milestones, both at the kinematical and dynamical level. First, using generalized power se- 
ries, purely kinematical definitions of these cosmological events are provided in terms of the 
behaviour of the scale factor a(t). The notion of a "scale- factor singularity" is defined, and 



143 



its relation to curvature singularities (polynomial and differential) is explored. Second, dy- 
namical information is extracted by using the Friedmann equations (without assuming even 
the existence of any equation of state) to place constraints on whether or not the classical 
energy conditions are satisfied at the cosmological milestones. We use these considerations 
to derive necessary and sufficient conditions for the existence of cosmological milestones such 
as bangs, bounces, crunches, rips, sudden singularities, and extremality events. Since the 
classification is extremely general, and modulo certain technical assumptions complete, the 
corresponding results are to a high degree model- independent: In particular, we provide 
a characterization of the class of bangs, crunches, and sudden singularities for which the 
dominant energy condition is satisfied. 

gr-qc/0508045; Accepted for publication in Classical and Quantum Gravity, 
celine . cattoenOmcs .vuw. ac .nz, matt . visserOmcs . vuw .ac.nz 

Note that the full version of this paper can be found in Class. Quant. Grav. 22:49 13- 
4930, 2005 or there is an e-Print Archive version: gr-qc/0508045. 
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Appendix C 

Gravastars must have anisotropic 
pressures 

Celine Cattoen, Tristan Faber, and Matt Visser 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

One of the very small number of serious alternatives to the usual concept of an astrophysical 
black hole is the "gravastar" model developed by Mazur and Mottola; and a related phase- 
transition model due to Laughlin et al. We consider a generalized class of similar models 
that exhibit continuous pressure — without the presence of infinitesimally thin shells. By 
considering the usual TOV equation for static solutions with negative central pressure, we 
find that gravastars cannot be perfect fluids — anisotropic pressures in the "crust" of a 
gravastar-like object are unavoidable. The anisotropic TOV equation can then be used to 
bound the pressure anisotropy. The transverse stresses that support a gravastar permit a 
higher compactness than is given by the Buchdahl-Bondi bound for perfect fluid stars. Fi- 
nally we comment on the qualitative features of the equation of state that gravastar material 
must have if it is to do the desired job of preventing horizon formation. 
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Appendix D 



Effective refractive index tensor for 
weak-field gravity 

Petarpa Boonserm, Celine Cattoen, Tristan Faber, 
Matt Visser, and Silke Weinfurtner 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

Gravitational lensing in a weak but otherwise arbitrary gravitational field can be described 
in terms of a 3 x 3 tensor, the "effective refractive index". If the sources generating the 
gravitational field all have small internal fluxes, stresses, and pressures, then this tensor is 
automatically isotropic and the "effective refractive index" is simply a scalar that can be 
determined in terms of a classic result involving the Newtonian gravitational potential. In 
contrast if anisotropic stresses are ever important then the gravitational field acts similarly 
to an anisotropic crystal. We derive simple formulae for the refractive index tensor, and 
indicate some situations in which this will be important. 
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